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We propose an improved scheme to construct many-body trial wave functions for fractional Chern 
insulators (FCI), using one-dimensional localized Wannier basis. The procedure borrows from the 
original scheme on a continuum cylinder, but is adapted to finite-size lattice systems with periodic 
boundaries. It fixes several issues of the continuum description that made the overlap with the exact 
ground states insignificant. The constructed lattice states are translationally invariant, and have 
the correct degeneracy as well as the correct relative and total momenta. Our prescription preserves 
the (possible) inversion symmetry of the lattice model, and is isotropic in the limit of flat Berry 
curvature. By relaxing the maximally localized hybrid Wannier orbital prescription, we can form an 
orthonormal basis of states which, upon gauge fixing, can be used in lieu of the Landau orbitals. We 
find that the exact ground states of several known FCI models at v — 1/3 filling are well captured 
by the lattice states constructed from the Laughlin wave function. The overlap is higher than 0.99 
in some models when the Hilbert space dimension is as large as 3 x 10 4 in each total momentum 
sector. 
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I. INTRODUCTION 

The Chern insulated is the first and simplest exam- 
ple of a topological insulator. It is defined by a non-zero 
Chern number of the occupied bands. It exhibits an inte- 
ger Hall conductance similar to the integer quantum Hall 
effect but at zero overall magnetic field. Recently, sev- 
eral groups added interactions to the one-body topologi- 
cal (Chern) insulator problem and reportecpHZl strongly- 
correlated phases similar to the Fractional Quantum Hall 
(FQH) effect. These fractional Chern insulators (FCI) 
have a partially filled band with non-zero Chern number 
and develop, similar to the FQH case, strongly correlated 
topological phases at specific values of v of the filling fac- 
tor. 

At filling v = 1/q (q = 3,5 for fermionPEE, q = 2,4 
for bosons^), a phase similar to the Laughlin state^ has 
been observed in numerical studies through several sig- 
natures. First, the system has g-fold quasi-degenerate 
ground states at filling v — 1/q, separated by a finite en- 
ergy gap from higher excitations. Second, twisted bound- 
ary conditions drive spectral flow within the ground state 
manifold, showing a Hall conductance a xy = e 2 /(qh). 
Third, the excitations from the topological ground state 
resemble the FQH quasiholesj^ identified by the count- 
ing^ of low- lyin g le vels in the gapped energy and entan- 
glement spectra™ 

Despite the similarity of the emergent features, the FCI 
and the FQH effects are hosted by substantially different 
systems. The FQH effects are typically observed in two- 
dimensional electron gas. This continuum system allows 
a holomorphic description, and the variational trial wave 
functions can be characterized by the asym ptotic behav- 
ior when particles approach each other DEI i n sharp con- 
trast, the Chern insulators are defined on a discrete lat- 
tice, and it s continuum limit is different from the FQH 
effect JEHU! This obstructs attempts to build model wave 



functions using the asymptotic behavior. 

An open question is thus to what extent and in what 
form the FCI can be described by the FQH physics. In 
this paper we provide a generic scheme to compare FQH 
and FCI wave functions at finite size. In particular, we 
address the fermionic FCI phase at filling v — 1/3 and 
clarify its relation to the Laughlin state. 

The first step in this direction was made by Qi in 
Ref. [15] The one-dimensional maximally localized Wan- 
nier states are plane waves localized in the direction per- 
pendicular to propagation. They form an alternative 
basis in a topological flat band, and they resemble the 
Landau orbitals in the lowest Landau level (LLL) in the 
continuum. The central, elegant idea of Ref. [15] is to ex- 
ploit this similarity and transcribe FQH wave functions 
written in the second-quantized form to the FCI system. 
The two-dimensional (2D) fractional topological insula- 
tor (FTI) system was similarly analyzed using decoupled 
FQH states in the Wannier basis. 

Upon closer inspection, the original proposal suffers 
from several issues that prevent direct application of this 
scheme to make contact with existing numerical stud- 
ies. The formalism in Ref. [15] was built upon a smooth 
gauge in the continuum. This poses technical challenges 
for numerical implementation on a finite-size lattice, and 
conceptually, the one-dimensional (ID) maximally local- 
ized Wannier states are orthogonal to each other only 
in the continuum limit. The non-orthogonality at finite 
size spoils the translational invariance of the constructed 
many-body states. More importantly, the maximally lo- 
calized Wannier orbital to LLL mapping in Ref. [15] is 
a mapping between the basis states of two separate sys- 
tems. The two sets of basis states have independent gauge 
freedoms. Merely introducing a ID relabeling of the Wan- 
nier states is not enough to properly fix the mapping from 
a Chern insulator to the LLL. A naive implementation of 
the formulas in Ref. [15] results in variational states very 
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low overlaps (lower than 0.04 for a system of 8 particles 
on a 6 x 4 lattice) with the exact diagonalization ground 
states of the two-orbital model discussed therein. 

In this paper, we build upon Ref. [15] and provide an 
improved prescription for constructing FCI wave func- 
tions using gauge-fixed (non-maximally) localized Wan- 
nier states, with large overlap with the exact diagonaliza- 
tion ground states. Our procedure is defined on a lattice 
of finite size with periodic boundary conditions in both 
directions. We do not need gauge smoothing, and we 
trade maximal localization for the orthogonality between 
the ID localized Wannier states at finite size. We explic- 
itly fix the phase choice of the Wannier states so that it 
best matches that of the Landau orbitals. Our prescrip- 
tion keeps both the center-of-mass translational symme- 
try and the inversion symmetry at the many-body level, 
and it implements the folding rul# that accounts for the 
total momenta of the degenerate FCI ground states on 
a torus. We also show that, in our prescription, using 
the Wannier bases localized in either x or y directions 
produces the same set of many-body states in the limit 
of flat Berry curvature. 

As a numerical test, we construct the Laughlin state 
for several FCI models at filling v = 1/3. In all cases, 
the lattice Laughlin states have very high overlap with 
the ground states, and their entanglement spectra share 
the same gapped structure. The overlap for the best 
FCI models found so far 7 reaches 0.99 for a system of 8 
particles on a 6 x 4 lattice, of which the Hilbert space 
dimension is as large as 3 x 10 4 in each total momentum 
sector. Our results provide the first direct, quantitative 
evidence that the fractionalized phase in Chern insulators 
at v = 1/3 can be well approximated by the Laughlin 
state. 

The paper is organized as follows. In Sec. [II] we discuss 
the many-body translational symmetries^ of the torus 
FQH states in the continuum. We construct a recom- 
bined set of degenerate FQH states that allow a direct 
mapping to the lattice. Sec. |III| summarizes with the 
basic properties of the ID localized Wannier states on 
a lattice. Except for a few subtleties, this discussion is 
similar to the one in Ref. [T5J 

In Sec. |IV| describes the Wannier construction of the 
many-body FQH states on a lattice. We elaborate upon 
the phase choice of the ID localized Wannier states. We 
show that the resulting many-body states are translation- 
ally invariant in both directions. We also demonstrate 
that the counting of the q-fo\d states in each momentum 
sector agrees with the counting rule proposed in Ref. |9j 
Eq. ( 70 ) provides an explicit and ready-to-use formula to 



rewrite any FQH wave function in a given Chern insula- 
tor model. 

Sec. |V| shows that our prescription preserves inversion 
symmetry on the many-body level when such symmetry 



exists at the single-particle level. Sec. [VT] presents the 
numerical verification that the ground states of various 
FCI models are well approximated by the Laughlin state 
constructed using our prescription. We also examine how 



close our phase choice for the Wannier basis is from the 
absolute ideal by a brute- force optimization. Sec. |VII 
examines the isotropy of the Wannier construction. We 
prove that the many-body lattice states constructed us- 
ing the Wannier bases localized in either direction are the 
same in the limit of flat Berry curvature. The isotropy 
in the presence of curvature fluctuations is tested numer- 



ically. Sec. VIII concludes the paper and discusses a few 
future directions. 



II. FQH TRANSLATIONAL SYMMETRIES 

It is instructive to review the many-body translational 
symmetries of the FQH system in the continuum. We 
begin by recapitulating the main results of the PRlPby 
Haldane. 

We study the problem of N e electrons moving on a 
twisted torus with a perpendicular magnetic field. The 
twisted torus can be represented by a parallelogram with 
opposite edges identified. We set up a Cartesian co- 
ordinate system (x,y), with orthonormal basis vectors 
(e x ,e y ). The tilde over x and y emphasizes the contin- 
uous nature of these variables. We consider the torus 
spanned by Li = L\e v and L2 = £2% (see Fig. [I]). 56 
Here, L\ and L2 are the length of the two fundamental 
cycles of the torus, and the unit vector e v is defined by 
e v = sin e x + cos e y: where the twist angle G (0, 7r) is 
the angle between the two fundamental cycles. The rect- 
angular torus corresponds to = tt/2. We introduce the 
normal unit vector e z = e x x e y , and define the reciprocal 
primitive vectors 

Gi = 2i\e x j{L\ sin 0), G2 = 2i\(e y — cot e x )/L,2. (1) 

They satisfy G a • L5 = 27r5 a b for a, b G {1,2}. 

The torus is pierced by a magnetic field in the — e z di- 
rection, B = V x A = Be z with B < 0. The Hamiltonian 
of the interacting electrons reads 



H 



[—ihVi - eA(Yi 
2m 



i<j 



We denote by e < t he charg e of the electron. The mag- 
netic length is Ib = y/h/{eB). The total number of fluxes 
penetrating the torus has to be an integer, given by 
LiL 2 sin<9 = 27r/|7V We define N = GCD(7V e , N^), 
where GCD stands for the Greatest Common Divisor. 
Then p = N e /N and q = N^/N are coprime. 

The magnetic translation operator of a single parti- 
cle is defined by T(a) = e ~* a * K /^, where K = —ihV — 
eA(r) + eB x r is the guiding center momentum of the 
z-th particle. We pick Landau gauge A(x,y) = Bxe yi 
and impose periodic (not twisted) boundary conditions 
T(L a ) = 1, a = 1,2. The wave functions (x,y\j) of the 
Ns single-particle states in the LLL are given by 
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<l>j(x,y) 



2tt(j + nA^)- 



L 2 



. 7r L\e 



-x 2 /(2l 2 B ) 



(3) 



r 



where the state index j is an integer defined modulo N^. 

We can decompose the translation operator of the i-th 
electron, T^(a), into a relative part and a center-of-mass 
part, Ti(a) = T reM (a)T cm (a/Ae), where 



Here and hereafter throughout the paper, we use short- 
hand notation for products and summations: 



M 



M-l 



M 



M-l 



n-n 



5>E- 



(5) 



i=0 



z=0 



We are particularly interested in the following translation 
operators and we define the shorthand notations 



Now we are in a position to establish a concrete rep- 
resentation of the many-body translation algebra. The 
basis states in the LLL are \{j}) = |j ,ji, • • • , 3N e -i) 
with an implied anti-symmetrization. The action of the 
many-body translation operators reads 9 



r reM (a)=T i (a)nr j (- 1 ^), T cm (a) = n^(a). (4) ^ 



|{J}) = I{J 
T?J{j}) = \{3 



q}), T?j{j})- 
■i}>, Tyj{j})- 



i2^j/N 



-i2n-£j/N4, 



\{j}), 



(9) 



Here the sum J^j runs over all the particles, and the 
state \{j + I}) is obtained from \{j}) by shifting the j 
quantum number of each particle by I. As expected, the 
relative and the center-of-mass translations are related 
by the periodic boundary conditions, T* el = [(T* m ) q ]^ 

r r y el = [(W]t. 



rpX 

1 rel,', 



T re i,»(pLi), 



rel,i 
cm 



T r el,i(pL 2 ), 

T cm (L 2 /A^). 



(6) 



Thanks to the periodic boundary condition, we can drop 
the particle index i from the relative translation opera- 
tors, since T* M = [(T^ and T^. = [(T^)^ do not 
depend on the particle index. 

The operators T? el and Tj^, and of course T^ el being a 
power of TJ^, commute with each other and the Hamilto- 
nian #P We use them to block-diagonalize the Hamilto- 
nian into momentum sectors labeled by two-dimensional 
(2D) wave numbers (k x , Ky} in an N x Brillouin zone, 
defined by the eigenvalues^ 



rel 



cm 



y/N* 



(7) 



The operator T£ m commutes with Tf el and H, but not 



with T C V 



rpy rpX _ JiX rjiy - 

cm cm cm cm 



-i27rp/q 



(8) 



Therefore, the many-body energy eigenstates can be 
grouped into g-fold center-of-mass multiplets; the q states 
in each multiplet can be transversed by successive appli- 
cations of T£ m and they share the same energy and the 
value of n x and n y mod TV. 

We set the cyclotron energy to infinity and focus on the 
states in the LLL. In this approximation, the many-body 
Hilbert space is spanned by the occupation-number basis 
states (Slater determinants) constructed from the LLL 
states and the projected interacting Hamiltonian is 
given in Appx. [Aj We use the curly braces {•} to denote 
a list of quantum numbers for the N e electrons. For suc- 
cinctness, we always omit the electron indices from {•} 
and related expressions. 



A. Recombination of the g-Fold States 

It has been suggested that the FCI with \C\ = 1 on 
an N x x N y lattice corresponds to a FQH system with 
fluxJV = N x N y W^ We now look for an alternative rep- 
resentation of the many-body translational symmetries 
with a Brillouin zone commensurate to the lattice sys- 
tem, in preparation for the analysis of the translational 



symmetry of the Wannier construction in Sec. IV C We 
emphasize that we are still working in the continuum 
and we are merely providing another representation of 
the center-of-mass translation algebra. The integers N x 
and N y should be understood as a factorization of at 
this stage. 

Following Ref. 9, we define the integers No x = 
GCD(N e ,N x ), N 0y = GCD(7V e , Ag, and p x = N e /N 0x , 
P y = N e /N 0y , q x N x /N 0x , q y N y /N 0y . Obviously 
p Xl q x are coprime, so are p y ,q y - Less obviously, q x q y 
divides g, 9 and thus p,q x are coprime, so are p,q y . We 
define the translation operators^ 



(T x 

V- 1 - ci 



\Q/Qx 



Ry — (^cm 



(10) 



They commute with each other and the Hamiltonian H. 
Therefore, we can use (S x ,R y ), instead of (T^Tj^), 
to block-diagonalize the Hamiltonian into momentum 
sectors. Since (S x ) qx = (T^)" 1 ", the new set of 2D 
wave numbers defined by (S x ,R y ) takes value from an 
(Nq x ) x (No x N y ) Brillouin zone [as opposed to the NxN^ 
Brillouin zone of Eq. 

Notice that S x operates within the g-fold multiplet. 
Successive applications of S x break the multiplet into 
non-overlapping orbits. The states in each orbit have 
the same R y eigenvalue due to [S x ,R y ] = 0. The orbit 
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structure can be revealed by the eigenvalues of X^ m , as 
they are distinct for the q states. Plugging Eq. (10) into 
Eq. (l8|), we have 



- L cm°x — ^ cm e 



i2np/q x 



rpX JD 

A cm Ix y 



D rpx A2irp/(q/q x ) 
1 h '!J ^ cm e 



(ii) 

(12) 



From Eq. (11), every orbit of S x has the same length q x . 



The total number of orbits is thus q/q x . From Eq. (12), 



for r G [0 .. q/q x ), (^fm) r brings any given state in the q- 
fold multiplet to q/q x distinct states discriminated by R y 
eigenvalues, and thus belonging to q/q x different orbits 
of S x . This covers each of the orbits exactly once. 

We can recombine the g-fold states to form simultane- 
ous eigenstates of S x and R y . We start the construction 
from an arbitrary energy eigenstate in the g-fold mul- 
tiplet diagonal in (T? el , T c ^ m ), labeled by k x G [0 .. TV) and 
K y G [0 .. Ny). The usual, relative y- momentum which 
labels the eigenvalue of T^ el can be easily obtained as n y 
mod N. For the Laughlin case at filling v — 1/3, is 
just one of the threefold degenerate states on the torus. 
The q states can be regrouped into the orbits of S X1 la- 
beled by r: 

{(S x ) m (T* m y\*) | m e [0 .. &)}, r € [0 .. q/q x ). (13) 

We can recombine the q x states in each orbit into q x 
eigenstates of S x . Within the multiplet, ince (S x ) Qx = 



, we have 



q x _ e -i27rq x (K x +sN e )/(Nq x ) 



(14) 



Note that the values of e ^(^+sN e )/(Nq x ) are distinct for 
s G [0 .. q x ), thanks to p, q x being coprime. We define 
the q x x (q/q x ) states 

-j q x 

|tf;s,r) = _^e <2wm ^ +a ^/^)(5 a! ) m (7Sn) r |*), 

(15) 

for s G [0 .. (fc) and r G [0 .. q/q x ). These states are 
orthonormal, and they are simultaneous eigenstates of 
S x and R y in the (Nq x ) x (No x N y ) Brillouin zone} 
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S x \V\s,r) =e 
R y \^;s,r) =e 



-i2Tr( Kx +sN e )/(Nq x )^. ^ r ^ 
-i2ir(K y +rN e )/(N 0x N y )^. s ^ r y 



(16) 



The q x x (q/q x ) states are related by center-of-mass 
translations, namely 



T^;s,r) = |¥; S ,r+l), 

T c ^ m |*;s,r) = e-^K+^V^^js-l,, 



(17) 



In terms of amplitudes, the second equation above can 
be written as [using Eq. ([9|] 

e- i2 *^ j ' N +({j}\9]s,r) 

= e -i2*( Ky +rN e )/N+({.y^. s _ ^ ( lg ) 



Incidentally, we note that s,r) is periodic under s —> 
s + q x , but it acquires a phase when r — >> r + q/q x [since 
(T* m )«/«* = ^ as defined in Eq. pob] 



B. Translational Symmetries in Amplitudes 

We now consider the manifestation of the translation 
symmetries in the amplitudes s, r). The action of 

S x and R y on the basis states \{j}) reads 



S x \{j}) = \{j+q/q x }), 
Ry\{j})=e- i2 ^ j/(No ^ ) \{j}). 

Therefore, from the first equation above we have 



(20) 



e- i2 ^+^)/(^)({j}|*;s,r) = ({j - q/q x }\V; s, r), 

(21) 

while from the second, we find that the amplitude 
s, r) vanishes unless 



j = *y + rN e mod N 0x N y . 



(22) 



The last two equations above summarize the information 
that we need from the FQH side to establish the trans- 
lational invariance of the FCI many-body wave functions 



to be constructed in Sec. IV on an N x x N y lattice. 



III. HYBRID LOCALIZED WANNIER STATES 

We now proceed to establish the hybrid Wannier basis 
in a Chern band on a lattice. The ultimate goal is to con- 
struct ID localized Wannier states which are plane waves 
in the second direction (hybrid). Such states mimic the 
Landau orbitals \j). We begin by reviewing the construc- 
tion and the properties of the ID maximally localized 
Wannier statesP^"^ 

Consider a 2D band insulator with orbitals per 
unit cell, indexed by a. We assume lattice translational 
symmetry and periodic boundary conditions. Denote 
the primitive translation vectors by bi and b2, with 
&z ' (bi x b2) > 0. Then the Bravais lattice is indexed 
by xhi + yb2, with (x, y) G Z 2 (be to differentiated from 
(x,y) G M 2 in the continuum.) We pick the principal re- 
gion to be (x,y) G [0 .. N x ) x [0 .. N y ). The momentum 
space is given by the reciprocal lattice k. We label points 
in the momentum space by wave numbers (k x ,k y ) G Z 2 , 
defined by k x = k • bi and k y = k • b2. The single- 
particle orbitals can be written as \x,y,a), or \k x ,k y ,a) 
in the momentum space. 

The energy eigenstates are Bloch waves. We focus on 
a single, isolated band \k x ,k y ), which will be fractionally 
filled in the many-body construction. For now we only 
focus on the one-band problem. Generalization to the 
2D FTI problem is straightforward as that problem de- 
couples in the Wannier basis. The wave function of the 
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Bloch band is u a (k x , k y ) = (k x , k yi a\k xi k y ). We assume 
periodic boundary condition in both k x and k y directions, 
\k x + l x N x , k y + l y Ny, n) = \k x , k y , n), for l x , l y e Z. We 
set the first Brillouin zone to [0 .. N x ) x [0 .. N y ). 

We denote the position of an orbital a relative to its 
unit cell coordinate by e^bi + e^,b2, where (e£, e y a ) Gl 2 
are the relative displacements of the orbitals within an 
unit cell. Taking into account the periodic boundary, we 
can define the (exponentiated) position operators 

* = E E I*' »' a)e-^(^)/^ y, a\, (23) 

and similarly y. Here the integers (x, y) are summed over 
the Bravais lattice points in the principal region of the 
torus, and a is summed over all A/5 orbitals. 
Using the projector to the occupied band 

P= \kx,ky)(k x ,ky\, (24) 

kx iky 

we can define the projected position operators X = PxP 
and y = PyP. In the momentum space, they take the 
form 

X = ^ ] |fca;, &2/)*Ae(Me) ky)(kx H~ lj 

^ ^ (25) 

3^= ^ ^ l^a;, ky)J\,y{k X ^ ky)(k X , ky + l|, 

where (k x ,k y ) is summed over the first Brillouin zone, 
and A Xj y(k x ,k y ) are the (exponentiated) Berry connec- 
tions given by 

*^-cc (Me 1 k"y) (k x , |x| /c^ H - 1, 
A y ( ^a: i^y) = { i k>y \ V \ ^x i ky ~~h 1 ) 

= Y, e ~ i2 ' eVjN "<^ K)^^, ky + 1). 

a 

(26) 

Notice that the Berry connections depend on the embed- 
ding of the lattice and its orbitals in real space, as does 
the Berry curvature.^ The link to the usual continuum 
definition of the Berry connections is explicitly demon- 
strated in Appx. [F] We emphasize that the definition of 
the exponentiated discrete Berry connections does not 
require a smooth gauge. 

A. Localization and Orthogonality 

Thanks to the translational invariance, we can view 
the 2D system as N y decoupled ID subsystems labeled 
by the momentum k y . The position operator X oper- 
ates within each subsystem and is diagonal in k y . At 
each value of k y , the eigenstates of this operator are the 



closest analogue to the states \x,k y ,a) that we can con- 
struct within the occupied bands; they are the maximally 
localized Wannier state for the ID subsystem EHUD The 
eigenvalue of X is the (exponentiated) Wannier center 
position in the bi direction. 

Ref. [15] proposed to build FQH wave functions using 
the maximally localized Wannier orbitals as one body ba- 
sis with a one-to-one mapping into the LLL orbitals. A 
major issue at finite size, however, is that the N x eigen- 
states of the projected position operator at each value 
of k y are not orthogonal due to the non-unitarity of the 
projected position operator at finite N x . This point is 
demonstrated in Appx. |B} 

Gram-Schmidt orthogonalization is not a suitable solu- 
tion, as it mixes the Wannier states localized in different 
unit cells, thereby spoiling the translational invariance. 
The proper resolution is to keep only the phase part of 
the exponentiated Berry connection and use the alterna- 
tive, unitary projected position operators to define the 
ID localized Wannier states. Specifically, we define the 
unitary connections^ 

A a {k x , ky) = A a {kxi k y )/\A a (k x , k y ) | , a = x,y, (27) 
and the unitary projected position operators 

X= ^ ] \k X j k y )A x (k x , ky)(k x + 1, k y \, (28) 

kx ^ky 

Y= ^ ^ \k X j k y )A y {k x , k y )(k x , k y -\- 1|. (^9) 

kx 1 ky 

Technically, the eigenstates of these operators are not 
maximally localized (they are almost so), but this is a 
small price to pay for the orthogonality. 

If there is more than one occupied band, such in the 
case of the quantum spin Hall (QSH) effect and 2D 
topological insulators, the prescription for obtaining or- 
thonormal orbitals is to perform a singular value decom- 
position (SVD) of the non-Abelian exponentiated Berry 
connection into VEW^, and use VW^ in place of the 
Berry connection to define the projected position opera- 
tors.^ Specifically, if the occupied bands are denoted by 
\k x ,k y ,n), we replace the matrix (k x , k y: m\k x + 1, k y , n) 
(of indices m, n) by its SVD when defining the projected 
position operator 

Thanks to the gauge periodicity in both k x and k y , we 
can define the unitary Wilson loops 

N X ^ 

W x {k y ) = J^J A x (k x , ky), W y {k x ) = J^J A y (k x , k y ). 

kx ky 

For later convenience, we introduce X x (k y ) and X y (k x ) 
defined by 

[Mky)] N ' = W X {ky), [\y(k X )] N » = W y (k X ). (30) 

We pick X x (k y ) to be the A^-th root with argument 
angle in (—2^/^,0], and X y (k x ) to be the A^-th root 
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with argument angle in (— 27r/X 2y , 0]. We can interpret 
X y {k x ) and X y (k x ) as the "average" phase (connection) 
of the Wilson loop evenly distributed to each bond along 
the loop. In general, for more than one occupied band, 
X x (k y ) and \ y {k x ) are obtained from the sets of eigenval- 
ues of the corresponding Wilson loop matrices. 

The eigenvalues and the eigenstates of X can be found 
easily, due to its decoupling into N y blocks labeled by k y . 
Each block X^ y is an N x x N x matrix in the Bloch basis. 
The N x eigenvalues of Xk y are given by the A^-th roots 

Of W X (ky)^ 



X ky \X,k y )=A X (X,k y )\X,ky), 



(31) 



where A x (X,k y ) = e~ i27rX / Nx \ x (k y ), indexed by X G 
[0 .. N x ). The eigenstate belonging to the eigenvalue 
A X (X, k y ) can be written 




Yl K x A x (k, k y ) 
Recall from Eq. ([5| that we use the shortand 



k x -, ky) . (32) 



k x 



| A x (a€, ky) — J^J A x (k,, ky). 



(33) 



Since the curly-braced prefactor in Eq. (|32| is unitary, 
these states inherit at finite size the orthonormality of 
\k x ,k y ) through the Fourier transform. These orthog- 
onal Wannier orbitals are not maximally localized for 
the finite-size lattice, but quickly become so as the num- 
ber of sites N x increases. They exhibit a gauge free- 
dom, $>y(X, k y ) G M, which is a phase to be specified in 
Sec. ITVBl 

The center of the ID localized Wannier state is given 

by 



(X, k y \x\X, ky) = A x (X,k y ) 



1 Nx 

y 1 \A x (k x , k Vj 



(34) 

where the bracketed factor is a positive real number. 
Therefore, we can interpret 



X 



X,k v 



arg[A,(X, ky)] =X~^ Big[W x (k y )] 



(35) 

as the center of the maximally localized Wannier func- 
tion, where we pick the branch of the argument function 
with arg(z) G (-2tt, 0] to make x X,ky € [X,X + 1). We 
place the X-th unit cell over [X, X + 1). As expected, in 
each unit cell there are N y Wannier centers, correspond- 
ing to the N y values of k y . 

For a flat-band Hamiltonian with no band dispersion, 
such as the single-particle part of the FCI Hamiltonians 



in Ref. 7J the Wannier states are actually single-particle 
energy eigenstates of the lowest, fractionally filled band. 
This qualifies the set of all the ID localized Wannier 
states {|X, k y )} as a single-particle basis for construct- 
ing the many-body trial states for the FCI, which lie en- 
tirely in the lowest band, in the same spirit as the LLL 
projected FQH wave functions. 



B. Gauge Freedom 

In the language of Ref. 18, the construction of eigen- 
states in Eq. ( [32] ) can be interpreted as transforming the 
Bloch states with an arbitrary phase factor to the "par- 
allel transport" (pt) gauge, in which 



I k x , k y / 



P t 



[K(ky)] 



Yl K x A x (k : k y ) 



I k x , ky ) , 



(36) 



and the ID localized Wannier states are just Fourier 
transform of the Bloch states \k x ,k y ) P f Despite the 
name, technically the phase of \k x ,k y ) pt contains not 
only the parallel transport Yl^y A x (k, k y ) from |0, k y ), it 
also contains a rotation to accommodate the Wilson loop 

W X (ky). 

Under the gauge transform 



\k x ^ky) — > e ir] ( x ' v ^\k x ,ky) , 
the Wilson line transforms by 

kx kx 

H A x (k, k y ) -»• e - ij?(0 ^ A x (k, k y ) e 



ir](k x ,k y ) 



(37) 



(38) 



and thus the Wannier state transforms by 

\X,k v )^e^°^\X,k v ). (39) 

The reason for the appearance of e lTi ^ kx ' kv) with only 
k x = is obvious: the parallel transport gauge [Eq. ( [36] )] 
depends only on the initial condition at k x = and the 
gauge- invariant Wilson loop W x (k y ). Had we chosen to 
start the parallel transport from k®, the Wannier state 
would acquire a phase e %r) ^ kx ' kv) upon the gauge trans- 
form. The full gauge freedom of the Wannier states is 
kept in the phase e l ^ v ^ x ' ky \ Due to the X-dependency, 
this is much larger than the gauge freedom in e^ '^). 

In Ref. 15, a similar gauge transformation formula was 
derived in a different, smooth gauge in the continuum 
limit that is not periodic in the k y direction. However, 
as detailed in Appx. |F3| we find in that case a gauge 
transformation formula that is fundamentally different 
from the one given in Ref. [151 due to the absence of 
periodic boundary condition in the k y direction in the 
setup of Ref. 15, a gauge transformation with a non- 
trivial winding may shift the ID localized Wannier state 
to another unit cell, changing the X index. 
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X A\ 



FIG. 1: Single-particle orbitals in the lowest Landau level 
(LLL) and the Chern band. In the left panel we show the 
Landau orbitals in the LLL on a torus. The two fundamen- 
tal cycles are marked by Li and L2, and the twist angle is 
labeled by 0. The Landau orbitals are plane waves in the e y 
direction that are localized in the direction perpendicular to 
propagation. In the right panel we show a ID Wannier state 
localized in the bi direction and in the lowest band of the 
Kagome lattice model. The size of the spheres depicts the 
weights of the Wannier state on each lattice site. 



Incidentally, we point out that any gauge fixing in the 
k y direction at k x ^ is futile, as the "parallel trans- 
port" construction would effectively override any gauge 
choice at k x 7^ 0. This is reflected by the independence 



of Eq. ([39]) from e ^'*^ with k x ^ 0. 



We emphasize that the phase factor e l ® y ( Xyky \ or the 
single-particle gauge for \k x ,k y ), do not affect the local- 
ization of the Wannier state, but they do affect essentially 
the similarity with the L LL. We will discuss the choice 
ofe i* y (x,k y ) in Sec.llVB 



C. Connection to the Lowest Landau Level 

Ultimately, we want to bridge the physics of a partially 
filled band with a Chern number C = ±1 and that of the 
spinless FQH effects in the LLL on a continuum torus. 
(Unless noted otherwise, in this paper we specialize to 
Chern number \C\ = 1.) 

As illustrated in Fig. [T] at the single-particle level, 
both the Wannier states \X,k y ) in a Chern band and the 
Landau-gauge orbitals \j) in the LLL are plane waves 
localized in the direction perpendicular to wave propa- 
gation. Building upon the correspondence between the 
center lines of the localized states, Ref. [15] proposed a 
linear relabeling of the Wannier states. The resulting ID 
label is analogous to the LLL state index j. We now 
examine this relabeling of the Wannier states in more 
details. As we will see, the original prescription lacks a 
few ingredients that are essential for a concrete numerical 
implementation of the scheme. 

First, we need to specify the torus geometry in the 
continuum. The number of fluxes in the corresponding 




FIG. 2: The principal Brillouin zone for C = +1. Plotted 
are the flows of the Wannier center position x X,ky — X — 
&rg[W x (k y )]/ (2n) [Eq. ( |35[ )] as a function of k y , color-coded 
according to the unit cell index X. The movement of the 
Wannier center in each unit cell is monotonic in the principal 
Brillouin zone Ck y +5 y G [0 .. N y ) (marked by the gray shade), 
but not monotonic in the Brillouin zone k y G [0 .. N y ). 



LLL problem is given by = N x N t 



Mm 



We want the 



continuum torus to take the same shape as the lattice 
system with periodic boundaries. The aspect ratio L\jL 2 
is set tcP 



Li 
L 2 



Ny\h 2 [ 



(40) 



The only remaining parameter is the twist angle be- 
tween Li and L2 (defined in Sec. [ll]). For brevity, we 
want to map both C = 1 and C = — 1 Chern bands to 
the LLL with a magnetic field in the — e z direction as 
studied in Sec. [TT[ For C = 1, we put Li and L 2 in the 
direction of bi and b 2 , respectively. For C = —1, we 
need to "flip" the fundamental parallelogram. This cor- 
responds to choosing the twist angle = tt — (bi,b2), 
where (bi,b 2 ) G (0,7r) is the angle between bi and b 2 . 



This choice of 6 is justified numerically in Sec. VI B 



We now make a detailed comparison between the LLL 
orbitals \j) and the Wannier orbitals \X,k y ). In the LLL, 
the center coordinate Lij/N^ of orbital \j) in the bi di- 
rection is a monotonically increasing function of j, and 
the j = orbital is centered along the line x = 0. In 
the Chern band, the unitary Wilson loop W x (k y ) is a 
pure phase; when k y changes, it winds around the unit 
circle at the origin of the complex plane. In the contin- 
uum limit, W x (k y ) has winding number C = ±1 when 
k y goes around a single Brillouin zone. Assume for the 
moment that the winding motion of W x (k y ) on the unit 
circle is unidirectional (clockwise when C = 1 and coun- 
terclockwise when C = — 1) when k y increase from to 
N y , i.e. there is no "zigzag" pattern of going back-and- 
forth anywhere. In this case, when the k y Brillouin zone 
boundary is properly chosen, aig[W x (k y )] G (— 27r, 0] de- 
pends monotonically on k y . 

The proper choice of the Brillouin zone is illustrated 
in Fig. |2j We start the k y Brillouin zone from the point 
where W x (k y ) is closest to 1 in the lower half of the com- 
plex plane. The specific prescription is the following. 
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We pick the branch of dJcg\W x (k y )\ that takes value in 
(— 27r, 0]. We define the shift S y G [0 .. N y ) as the cardi- 
nality of the set 

{ ky G [0 Ny) | 3Ig[W X (ky)] > ^[W^O)] }. (41) 

Then, for Ck y + S y G [0 .. TV^), arg[W x (&^)] depends 
monotonically on k y ; Ck y + 5 y labels the Wannier cen- 
ters in an unit cell sequentially in ascending order of the 
center position x X,ky • We refer to the set of k y given by 



Cky+Sy G [0 .. Ny) 



(42) 



as the "principal Brillouin zone" (pBZ). We introduce 
the ID label of Wannier states 



jX,ky = X M + Ck V 



(43) 



with Ck v 



S y e [0 



Ny). We need to shift k y back 



to the pBZ before performing the mapping. Thanks to 
the explicit choice of the domain for k y , the mapping 
(X, k y ) — >• j is invertible. In the presence of "zigzag" 
patterns in W x (k y ), we stick to the above j x > k v formula, 
and take S y as a variational parameter. This simple de- 
pendence on k y is necessary for the translational invari- 
ance of the FQH-analog wave functions, as we will see in 
Sec. IV C| In the rest of this paper we will focus on the 
case without "zigzag" patterns. 

The phase factor e l ^ y ^ x,ky " ) is subject to periodicity 
constraints. Since in the LLL \j) = Ij+TV^), the Wannier 
states \X,k y ) and \X + N x , k y ) are mapped to the same 
LLL state. This mandates the gauge periodicity in X, 
i.e. \X, k y ) = \X + N x , k y ). Eq. ([32]) then leads to 



e i$y(X,k y ) = e i$>y{X + N x ,k y ) _ 



(44) 



The periodic gauge of the Bloch state in the k y direction, 
i.e. \k x ,k y ) = \k x , k y + N y ), requires that 



e i®y(X,ky) = e i<S>y{X,ky+Ny) _ 



(45) 



It should be noted that the analogy between the Wan- 
nier states and the LLL states is not exact even in the 
continuum limit, due to the inevitable fluctuations in the 
Berry curvature. For example, if the Berry curvature is 
not uniform in the y direction, the Wannier centers x x,ky 
are not distributed evenly over each unit cell, in contrast 
to the uniform distribution of the LLL orbitals. This is 
in line with the comparison between the LLL and the 
Chern band made in Refs. I20I2T1 



IV. WANNIER CONSTRUCTION OF FQH 
STATES ON A LATTICE 

We now turn to the construction of model wave func- 
tions in a fractional Chern insulator defined on an N x x 
N y lattice. We consider finding model wave functions 



that approximate the ground states of generic density- 
density interaction with translational invariance, namely 



Via 



r-h5,/x> 



(46) 



Here 6 is the unit cell displacements between neigh- 
boring orbitals \i and v. The sum (/i, z/, S) is over all 
distinct (up to lattice translation) orbital pairs within 
some interaction range. The Bravais lattice coordinate 
r = (x, y) G Z 2 is summed over the principal region of 
the torus. We refer the reader to Ref. 7 for a series of 
examples with concrete choices of (/i, S). Longer range 
interactions can also be implemented, but our experience 
is that they diminish the strength of the FQH-like state. 

In the flat-band limit of the single band approxima- 
tion^, the interacting Hamiltonian is just the projected 
density-density interaction PVi^P: 



lat 



L v v v „-i(k 1 -k 4 )-(5 



^ {fj,,u,S} kik 2 k 3 k 4 
X ^( k l)<( k 2)^( k 3)^( k 4)^k 1 ^k 2 ^ k3 ^ k4 . (47) 



Here k n = (k Xin , k y ^ n ), n = 1, 2, 3, 4, are summed over the 
first Brillouin zone, and the primed Kronecker-c) allows 
umklapp processes ki + k2 = k3 + k4 mod (N x , N y ). 

Based on the physics of the FQH, we expect that the 
Hamiltonian in Eq. (47) could have a topological ground 
state at filling v = 1 /qthat resembles the Laughlin state. 
The Laughlin state expansion in non-interacting many- 
body states is known and hence all we need for a model 
FCI state is an appropriate map between the FCI single- 
particle orbitals and those of the LLL. The relevance of 
these trial wave functions to the actual FCI ground state 
is demonstrated in Sec. EH We now examine the details 
of this construction. 

For a flat-band single-particle lattice Hamiltonian, the 
many-body Hilbert space is spanned by the Slater deter- 
minant states \{X, k y }). The ID index j x ^ k v defined in 
Eq. (43) provides a formal mapping between the many- 



body basis states |{X, k y }) on the lattice and \ {j}) in the 
continuum. The lattice analogue of |^; 5, r) states can be 
constructed: 



|*;s,r)i at = 



E \{X,ky}){{j XM 

{X,ky} 



}|*;*,r), (48) 



where ({j Xjky }\ i &; s,r) are the FQH amplitudes of the 
states defined in Eq. (15), and {X, k y } is summed over 
all A/" e -particle configurations in the Wannier basis. 

Naively, this seems to be the end of the story. How- 
ever, the state s,r)i a t defined above is not covariant 
under a single-particle gauge transform on \X,k y ), as the 
continuum states \j) do not transform accordingly. We 
have to fix the phase e l ^ y ^ x ' ky ^ of the Wannier states 
|X, k y ) in conformity with the phase of the \j) states. 
This is not surprising: following Ref. [15j up to now we 
have only established a mapping between the state labels 
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(X, k y ) and j, rather than a mapping between the actual 
states. In the following, we seek the guidelines for choos- 
ing e l ® y ( x > ky } based on the similarity with the LLL. If 
we do not properly fix the e l ^ y ^ x,ky ^ phase, the overlaps 
with the exact ground states can be tiny. 



A. Connection between Wannier States 

For the LLL orbitals \j) defined in Eq. ([3]), we find in 
Appx. [A] that the adjacent orbitals satisfy 



(j|e- iG2 - r |j + l) = e-l G2 l 2 ^/ 4 G 



(49) 



Here, G2 = 27r(e y — cot 6 e x ) j X 2 is the reciprocal lattice 
vector defined in Eq. in Sec. [TTJ and the position op- 
erator r takes values in the continuum. We can interpret 
e -%G 2 -r ag a translation by — G2 in the momentum space. 

This condition spells out the gauge choice of the LLL 
orbitals. In the continuum limit, it gives the paral- 
lel transport between the \j) and \j + 1) orbitals. We 
need to find a consistent gauge for the Wannier states. 
Recall that Y implements translation in the momen- 
tum space along k y , which is the direction reciprocal 
to b 2 . Therefore, on a lattice, the quantity analogous 
to (j\e~ lG2 ' r \j + 1) is the connection between adjacent 
Wannier states 



(X,k y \Y\X',k' y ). 



(50) 



Here, in the same spirit as the orthogonality fix for 
\X,k y ), we use the unitary projected position operator 



F, defined in Eq. ([29]), and (X' \k' y ) is the label of the 
Wannier state next to \X,k y ), defined by 



aX'X, 



C. 



(51) 



Plugging in the definition of j f ,ky in Eq. (|43j), we find 
k f y = k y + 1 mod N y . If going from k y to k y + 1 crosses 
the boundary of the principal Brillouin zone [Eq. ([42])] , 
the Wannier centers (X,k y ) and (X',k y ) are located in 
different unit cells (X' = X + C), otherwise they are in 
the same unit cell (X f = X). Notice that for \C\ > 1, 
Eq. (51 ) breaks the N x N y Wannier states into \C\ groups, 
and it implies a \C\ -layer FQH analogy.^ 

We ^can calculate the Wannier connection 
(X, k y \Y\X', ky) by expanding the Wannier states 

in Bloch basis using Eq. (32). Since the operator Y is 



diagonal in k x , the two Bloch basis expansions is reduced 
into a single sum over k x . Consider a generic term k x in 
the sum. The A x factors from the Wannier states and 
the A y factor from Y can be collected into a Wilson line 
going along (0,k y + l)(k x ,k y + l)(k x ,k y )(0,k y ), which 
equals A y (0,k y ) times the Wilson loop around the path 
shown in Fig. [3] The X x factors from the Wannier states 
are [X x (k y + 1) / X x (k y )) kx . If the two Wannier states are 
not in the same unit cell, they must be in adjacent unit 
cells X' = X + C [due to Eq. (51 )], and we have another 
factor e ^k x c/N^ 



k +1 

y 



BZ 



FIG. 3: The loop (0, fc„)(0, k y + l)(k x , k y + l)(k x , k y )(0, k y ) 
in the N x x N y lattice Brillouin zone. The Wilson loop around 
this path is defined in Eq. (54) as W^(k x ,k y ). Notice that 
at k x = N x , the connections at the two vertical edges cancel 
each other due to periodic boundary, leading to W^(N X , k y ) — 

W X (ky)/W X (ky + l). 



Putting all these together, for (X f ,k f y ) defined in 
Eq. (51), we have the expression 



(X,ky\Y\X'X) 



J*y(X',k') 



^5-^(0 A) (52) 



where the gauge-invariant quantity U y {k y ) is defined as 



1 ^ W_ 



W_(k x ,k y ) 



N X ^W.(k X ,ky) 

Here, W^(k x , k y ) is the unitary Wilson loop in Fig. [3j 

Tll^ A X (K,,ky) Ay(k Xj ky) 



(53) 



W_(k x ,k y ) 



Yl k K *A x ( K ,k v + l) A y (0,k y ) 



(54) 



while W^(k x , k y ) — [/^ x (k y )] kx is given by /J JX (k y ) defined 
over the principal Brillouin zone [pBZ, the set of n satis- 
fying Ck + J„ G [0 .. N v )], 



l^x(ky) 



^X\ky + 

J2ttC/N x 



^x{ky) 

X X (ky + 1) 



if k y + 1 e pBZ, 



otherwise. 



(55) 



We emphasize that the above definition requires first 
shifting k y back to the pBZ. The two cases in Eq. (55) 
correspond to whether or not going from k y to k y + 1 
crosses the boundary of the pBZ. 

We now try to understand fi x (k y ) and W^{k xi k y ) 
physically. Without loss of generality, we consider the 
case of C = +1. Since [X x (k y )Y^ = W x (k y ) [Eq. Q], 
we have for both cases of Eq. (|55|) , 



\p>x(k y )] 



N x 



X X (ky + l)_ 



W_(N x ,k y ). (56) 



Now we need to take the iV^-th root of this equation. Spe- 
cial care needs to be taken with the branch choice. Recall 
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a > &rg[\ x (ky + 1)] &rg[\ x (ky)} 



2tt 

' N x 



b) 



J arg[A x (^); 
aig[X x (k y )] &ig[X x (k y + 1)] 



2tt^ 



arg[A x (ft)] 



FIG. 4: Flow of arg[\ x (K,)] in the interval (-2ir/N x , 0]. The 
solid dots represents the X y values of arg[A x (ft)] for the X y 
Wannier centers in each unit cell. Here we show the case 
of C = +1: arg[A x (ft)] is a monotonically decreasing func- 
tion of k in the pBZ. Panel a) When k y and k y + 1 can be 
put in the principal Brillouin zone simultaneously, we have 
arg[X x (ky)/X x (k y + 1)] G (0,2n/N x ). Panel b) When go- 
ing from k y to k y + 1 crosses the boundary of the princi- 
pal Brillouin zone, we have <ng[e t2ir/Nx X x (k y )/X x (k y + 1)] G 
(0,2tt/N x ). Summarizing both cases, we have aig[/j JX (k y )] G 
(0, 2tv/N x ), for /^(fcy) defined in Eq. ([55]). 



from Sec. Ill A that arg[A x (ft)] lies in (—2tt/N x , 0] by def- 
inition. Since we assume the absence of "zigzag" patterns 
in the winding of W x (k), the motion of arg[A x (ft)] in the 
interval (—27r/N x ,0] must be monotonic when ft G pBZ, 
and arg[A x (ft)] jumps from the left boundary of the in- 
terval to the right boundary when k crosses the bound- 
ary of the pBZ. As illustrated in Fig. |4j for all values 
of k y , we have arg[/i £C (fc 2/ )] G (0,2tt/N x ). On the other 
hand, since the phase angle of the Wilson loop is given by 
the curvature enclosed in the loop, the argument angle 
arg[W_(iV £C , fey)] G [0, 2ir) is the total curvature enclosed 
in the N x plaquettes in the row between k y and k y + 1. 
Then, the N x -th root of Eq. (56) gives 



axg[n x (k y 



-arg[WL(JV x ,fc y )], 



(57) 



Therefore the argument angle arg[iy_(fe iC , k y )\ = 
k x a 1 Tg[ h L x (ky)] is k x /N x times the total curvature enclosed 
in the N x plaquettes in the row between k y and k y + I. 
When the curvature is constant across the N x plaquettes 
at k y , we have W^(k Xl k y ) = W^(k Xl k y ). The phase of 
each term in U y (k y ) [Eq. ( [53] )] can hence be interpreted 
as a measure of curvature fluctuations in the k x direction. 



B. Phase Fixing: Explicit Prescription 

Compared with {j\e~ iG2 '^\j + 1) G R+ in the LLL 
[Eq. p9])], a major difference in the Chern band is 
that the connection (X, k y \Y\X' , k') between the Wan- 
nier states cannot be gauge fixed to a real number. 
To see this, we take the product of Eq. ( [52] ) over all 
(X, k y ). Thanks to the periodicity of e^ (x '^[Eqs. (|44| 



and (45)], the product takes a simple form: 



N X Ny 

HH(X,k v \Y\X',k' y ) 

X k v 



N v 



Wy(0)l[U y (k y ) 



, (58) 



Without further constraints from symmetry, the expres- 
sion on the right hand side is a complex, rather than 
a positive real number. Even if we ignore the curva- 
ture fluctuations, in general we still cannot make all 
(X,k y \Y\X',k') G R+. This is fundamentally different 
from the LLL [Eq. ([49])]. 

Hence, we relax this condition and look for e l ^ y ^ x,ky ^ 
that makes the phase of (X, k y \Y\X' , k' y ) independent 
from X and fe^P^This phase is given by the A^-th root 
of the right hand side of Eq. ( [58] ) . Define the phase 
U y (k y ) = Uy(ky)/\Uy(k y )\. Similar to the phase X y (k x ) 
for the connection A y (k x ,k y ) [Eq. ([30])], we can define a 
phase uo y with argument angle in (—7r/N y ,7r/N y ] by 



N v 



{Uy) Ny 



l[U y (K). 



(59) 



Then, to have the phase of (X, k y \Y\X' , k' ) independent 



from X and k y , we can use 



(X,ky\Y\X',k' y ) 



\y(0)UJy\U y (ky)\. 



(60) 



Here, we have picked a specific branch of the N^-th root 
of (the phase of) Eq. (58). We could have picked an- 



other gauge by putting on the right hand side of Eq. ( 60 ) 
an extra factor e l27rw / N ^^ with w G (0 .. N^). However, 
this makes no difference to the resulting set of many-body 
states {|^; s, r)i at } [Eq. (48)]. The reason is the following. 
Without loss of generality, we consider the effect of this 
extra factor e l27VW / N ^ in the case with C > 0. When we 



gauge fix the Wannier states one by one, using Eq. (60) 
but with an extra e l27rw / N ^ between each pair of adjacent 
Wannier states, the Wannier state \X,k y ) receives an ex- 
tra factor e i27rwj ' y / N *. Therefore, the many-body state 



in Eq. (48) becomes 



{X,ky} 

e ; 2 ™ E i x ' k - (yx,k v }|$ ; S) r ), (6i) 



Plugging in Eq. ([18]), we find 



2-K W ( Ky +rN,)/N 4> ^. s + w ^^ 



lat- 



(62) 



The ^-dependence is reduced to an overall phase factor, 
and a shift in the s index, which can be absorbed by a 
reshuffling of s. Therefore, we can safely set w to zero, 
and use Eq. (60) to fix the phase of the Wannier states. 



Plugging this prescription into Eq. rt52|, we find 



e i<s> y (x f ,k y )-i$ y (x,k y ) = V°) 



Ay(0, ky) Uy(ky) 



(63) 
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Notices that the quantities on the right hand side are all 
unitary. Together with the initial condition e l ^y( ^) = i ? 
Eq. (63) recursively specifies the choice of e l ^ y ^ x,ky ^ for 



all the N(k states. Thanks to Eqs. (30) and (59), this 



choice of e i ® y ( x,ky ^ does not depend on X, and we can 
drop the X argument and write e l ^ y ^ ky ^ instead. We note 
that any choice of e i ^ y ^ ky ^ without an X-dependence that 
is periodic in k y can be achieved by simply modifying the 
gauge of the single-particle Bloch states along k x = 0. 
This changes the phase of the Wannier states by e l7] ( 0,ky ^ 
as shown in Eq. (39). 



Translational Invar iance 



specialize to the phase choice in Eq. (63) and consider 



the constraint from translational invariance on a generic 

Define the 7V e -particle center-of-mass translation oper- 
ators on the lattice 



rriX 

lat, cm 



{x,y,a} 



\{x + l,y,a})({x,y,a}\, 



7?at,cm= E \i X >V 



l,a})({x,y,a}\. 



(64) 



We now examine the translational symmetry of the 
many-body states \^>;s,r). For the moment we do not 



J 



The action on the wave function of the many-body lattice state \^>; s, r)i at defined in Eq. (48) is found to be 



<{X,fcJ|T£ tiC J*; S ,r) lat = e* 
({X, k v }\TV\V; s, r)ut = e~ i2 ^ ({X, ky}\9; s, 



r)i a t, 



(65) 
(66) 



where means summation over all the N e phases of different quantum numbers {X, k y } of the N e electron wave 
function. The translational invariance in the y direction is already apparent: Eq. (22) dictates that the non- vanishing 
components of ({X, s, r)i at have the same value of ^ k y = C[n y + (r — 5y)N e ] mod N y . This is the total 

momentum K y . The prefactor in Eq. (66) is the same for all the non- vanishing components, i.e. the wave function 



is translationally invariant in the y direction. The situation in the x direction is more involved. Plugging Eqs. (21) 
and (43) into Eq. (65), we have 



({X,fcJ|T 1 ^ cm |*; S ,r) lat =e- 



-i2ir(K x +sN e )/N x p iJ2l<S>y(X,k y )-<S> v (X-l,k y )] 



({X, k y }\^;s,r) iat . 



(67) 



In general, the exponential prefactor is different for 
each component (Slater determinants of different sets 
{X, k y }), spoiling the translational invariance. To restore 
this symmetry, e ' l ^y( x ^ k y)-' l ^y( x - 1 ^ k y) nas to be indepen- 
dent from X and k y . Translational invariance is part of 
the main reason why we also have asked for the phase of 
(X, k y \Y\X' , k f y ) to be independent from X and k y in the 
previous section. The periodic boundary condition on 
e i<5>y(x,k y ) then guarantees that the exponential prefac- 
tor is an A^-th root of unity, and thus the state recovers 
translational invariance in the x direction. 



D. Total Momenta and the Folding Picture 



We now consider the phase choice given in Sec. IV B 
Plugging Eq. (63) into Eq. (|67[), we have 



({X,fcJ|Tf atiCm |*; S ,r) lat 

= e- i2 ^^ + ^)/^({X,fcJ|vl>; S ,r) lat) (68) 

where the relative momentum n x of the FQH state is 
defined by the eigenvalue of the relative translation oper- 



ator T r gj in Eq. Hence the FQH-analogue many-body 
wave functions ({X, k y } \ ^; s, r) i at constructed through 
the Wannier functions using the above prescription are 
indeed translationally invariant. The total momentum of 
the 7V e particles in state s,r)i at is given by the wave 
numbers 



Kx — k x + sN e mod N x , 

K y = C[n y + (r - S y )N e ] mod N y . 



(69) 



Ref.|9] obtained the counting for FCI model states from 
the corresponding counting for FQH model states by 
first folding the N x N relative Brillouin zone down to 
No x x No y , and then unfolding to the N x x N y lattice Bril- 
louin zone. We now show that our procedure precisely 
reproduces this folding picture. 

First consider the case C = +1. Recall that 
GCB(N e ,N x ) = N 0x , N x = q x N 0x , GCD(N e ,N v ) = 
No y , Ny — q y No y . For s G [0 .. q x ), the q x values of 
K x are all distinct and can be written as k x + t x No x 
mod N x with t x G [0 .. q x ). This clearly implements 
the folding rule in the x direction. In the y direction, 
for r G [0 .. q/q x ), the q/q x values of K y form nothing 
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but a q/(q x q y )-fo\d replica of the q y values of K y +t y No x 
mod N y with t y G [0 .. q y ). This corresponds to the fold- 
ing rule in the y direction, producing q/(q x q y ) states in 
each momentum sector. The shift parameter S y can be 
absorbed into r and thus only reshuffles the order of K y 
values. For the case of C = —1, the minus sign from C in 
K y [Eq. ( [69] )] can be absorbed as a reshuffling of r. This 
is explained in details later, in Eq. (79). The results are 

thus in full agreement with Ref. [9j 

In Ref. [23l the generalized Pauli principl e 1 11 * 12 * was in- 
voked through the Wannier mapping^ to determine the 
total number of Moore- Read FCI quasihole states. Using 



our updated formalism, the number of quasihole states in 
each momentum sector could be found as well, and it is 
in agreement with the earlier results obtained in Ref. [23l 



E. Many-Body Amplitudes in the Bloch Basis 

We are now in a position to give the final formula 
for the amplitudes of the many-body lattice states con- 
structed from the FQH states s, r), 



({^,^}|^;s,r)i at = Yl 



{ujyfy [X y m ky [Ax(fcy)] 



Yt v U y {K) Ut y 4,(0,*) 



^M!! 1 _JL^ y e -«* E k x X/N x {{j X,k ym r y 



(70) 



Here the product outside the curly braces is over the (k x , k y ) configurations of the N e particles, and ({j x,ky }\^/] s,r) 
are the many-body amplitudes of the recombined FQH state^R _ ^ _ _ 



The above formula is the central result of this paper. This 
prescription is not limited to a specific model wave func- 
tion. It applies to any FQH state that can be expressed 
in a second-quantized basis; a real-space wave function 
is not actually needed. We note that the absolute value 
of each component in the Bloch basis does not depend 
on the gauge choice of the Wannier states, thanks to the 
X-independence of e l ® y ( x > ky \ 



V. INVERSION SYMMETRY 

Although not a vital element of a topological phase, the 
inversion symmetry is respected by FQH model wave- 
functions. Many Chern insulator models are inversion 
symmetric as well. For such models, the acceptable 
many-body wave functions for a featureless liquid state 
must respect inversion symmetry. To this end, in the 
following we specialize to inversion-symmetric Chern in- 
sulators and show that our construction always preserves 
inversion symmetry on the many-body level. 



A. Inversion Symmetry of the FQH States 

First we briefly discuss the FQH inversion symmetry in 
the continuum. Inversion is implemented by the unitary 
operator V: 



V\x,y) = \-x,-y). 



(71) 



Evidently V 2 = 1. Using the definition of (x,y\j) in 
Eq. (||, we can show that the annihilation operator \£j of 



state \ j) satisfies V^jV = ^-j, and thus the occupation- 
number basis states satisfy 



nm = K-j}>- 



(72) 



It is straight-forward to show that the FQH Hamiltonian 
H commutes with V, and that inversion flips the direc- 
tion of translation operators, 



VTV T\ 



(73) 



where T can be any of single-particle or many-body 
translation operator, such as X^, S x , or R y . 

Inversion does not commute with the relative trans- 
lation operators. Therefore, in general the operator V 
mixes different center-of-mass multiplets. However, the 
ground state energy level in an Abelian topological phase 
is populated exclusively by the g-fold degenerate ground 
state, and the operator V must operate within the g-fold 
center-of-mass multiplet, since [H, V] = 0. 

Consider the action of V on the state that we pick 
from the g-fold ground states. This state has total mo- 
mentum (n Xl n y ). By the above argument, n x must be in- 
version symmetric, while n y changes to — n y under inver- 
sion. Inversion on can be compensated by a center- 
of-mass translation. We can find a unique t G [0 .. q) 
such that 



tN P 



-2k v mod 



(74) 



Then, \^) must be an eigenstate of V(T^ m )K This 
composite operator satisfies [P(T^ m Y] 2 = 1 thanks to 
Eq. (73), and thus the eigenvalue has to be ±1: 



P|*) = C*(T ( f m ) t |*). 



(75) 
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We move on to discuss the recombined states s, r). 
Plugging Eq. (75) into Eq. (15), we have 



7>|tf;s,r) 



1 qx 

(76) 



/Qx 

=C*I^M, f), 

Here f is defined by 

f = t — r. (77) 

and s E [0 .. ^) by 

- (k x + siV e ) = k x + sN e mod 7V^. (78) 

It is easy to show that there exists a unique value of s G 
[0 .. q x ) satisfying the above constraint. From r = t — r, 
we can show that 

- (ay + r7V e ) = + f7V e mod N 0x N y . (79) 

In components, Eq. ( |76| ) reads 

({-j}|*; S ,r) = C*({i}|*;s,f). (80) 

B. Inversion Operator on the Lattice 

We now extend the inversion operator V to the lat- 
tice. As shown in Appx. |C 1[ the exponentiated position 
operators transforms under inversion by 



VxV = x 1 ", 



VyP = yl 



(81) 



Since the single-particle Hamiltonian is inversion sym- 
metric, the unitary operator V must take the Bloch state 
\k xi k y ) to |— k x , — k y ), up to a /c-dependent phase factor 



\~k X ,-ky) = e % ^^V\k X ,ky). 

Considering that V 2 = 1, we have 



(82) 



(83) 



These phase factors constitute the inversion sewing ma- 
trix; the explicit expressions can be found in Appx. |C 1| 
Since we work in periodic gauge, the phase factor must 
be periodic as well, 



(J-^kx+mNx ,ky+nN y — ^^k x ,k y 



(84) 



for m, n G Z. If both N x and N y are even, there are 
four points (k x ,k y ) in the Brillouin zone that are in- 
version symmetric (up to a reciprocal vector), namely 
(0,0), (JVa/2,0), (0,^/2), and (N x /2,N y /2). Thanks 
to Eqs. @ and ([MJ, we have e^ k ^ k v = ±1 if (k x ,k y ) 
is inversion symmetric, and e l ^ kx > k y at these points are 
just the inversion eigenvalues of the lowest band. The 
periodic boundary condition is crucial in relating e l ^ kx ' k y 
to inversion eigenvalue if (k x ,k y ) ^ (0,0). 



C. Inversion Symmetry of the Connections 

We now look for the implication of the inversion sym- 
metry in the Berry connection and the Wilson loop. The 
Berry connection satisfies 

A X (k X ,ky) = (k X ,ky\V 2 XV 2 \k X + l,ky) (85) 

= e^k x ,ky-^ kx+ i, ky (_ kx ^ -k y \x ] \-k x - 1, -ky) 



i£k x ,ky —i£k x + l,k 



y[A x (-k x -l,-k y )Y 



or for the unitary connections, 

A x {k xi k y )A x (—k x — 1, —ky) - 
A y {k x , ky)A y ( K —k Xl —ky — 1) : 



(A£k x ,k y i£k x + l,k y 

P ^k x ,ky —i£k x ,ky+l 



(86) 
(87) 



These are the lattice versions of the transformation of the 
Berry connections in the continuum [Eq. (73) in Ref. l24l] 
We now derive the lattice version of some well-known 
inversion properties of the Wilson loop. 24 Due to the pe- 
riodicity of the phase factors e z ^ kx > k y , the unitary discrete 
Wilson loop satisfies 



W X (ky)W X (-ky) = Wy(k x )Wy(-k x 



1. 



Therefore, at inversion symmetric momentum (k x ,k y ), 
we have 



W X (ky)=±l, 



Wy(k x ) = 



(89) 



At inversion symmetric k y , the Wannier center must be 
either on the Bravais lattice, or at a mid-bond point of 
the Bravais lattice. 

If N x is even, we can go farther and express W x (k y ) as 
a product of the two inversion eigenvalues at k y : 



.25 



N x /2 

W^xiky) ~ A X (k,, ky)A x ( Kj \,ky) 

— e l ^0,k y — z ^N x /2,k y ^ 

Similarly, if N y is even, we can show that 

WyCk X ) 



i£k x ,0 ^k 



x,N y /2 



(90) 
(91) 

(92) 



For the moment, we consider the continuum limit 
with large and even N x ,N y . Ref. l24l showed that the 
product of the four inversion eigenvalues at (k X: k y ) = 
(0,0), (A^/2,0), (0,^/2), and (N x /2,N y /2) is equal 
to ( — l) c = —1 for \C\ = 1. Therefore, in addition to 
Eq. (l88|), we now have 



W x (0)W x (N y /2) = Wy(0)Wy(N x /2) 



-1. 



(93) 



The case where N x or N y is odd needs special treatment, 
as we will see soon. 

These properties of the Wilson loop strongly constrain 
the shift parameter 5 y for the Wannier-LLL mapping in- 
troduced in Sec. |III C| Recall that we pick the argument 
angle arg[W x (k y )] G (-2tt,0]. The shift S y G [0 .. N y ) is 
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a) WJ0)=1 




S y =-(Ny-l)/2 



FIG. 5: The shift parameter 5 y fixed by inversion for C = +1. The three panels correspond to the three cases in Eq. (94). 
Compared with Fig. [2] here we highlight the discrete nature of the finite-size system. The wave number k y takes N y discrete 
values in each Brillouin zone. The solid dots represents the N y Wannier centers in each unit cell, colored coded by the unit cell 
index X. The gray shade shows the principal Brillouin zone Ck y + 5 y G [0 .. N y ). 



defined as the number of W x (k y ) with diig[W x (k y )] > 
arg[W x (0)]. We have shown that the Wilson loops 
W x (k y ) come in complex-conjugate pairs at momenta 
coupled by inversion, and in particular, (0) = ±1. 
As illustrated in Fig. [5j with inversion, the shift S y is 
fixed to 



Considering that we place the X-th unit cell over the 
interval [X, X + 1), this rule of transformation is not 
surprising: for example, inversion takes an on-site Wan- 
nier state centered at x = 2 + to a state centered at 
x = —2 + 0, while it takes a Wannier state centered at 
x = 2 + 0.3 to x = -3 + 0.7, to be in [X, X + 1). 



ifW^O) 
N y /2 ]£W X (0) 
(N y -l)/2 ]fW x (0) 



A key quantity in our gauge choice Eq. (63) for the 
Wannier states is \ y (0)ujy. We now examine the con- 
straint from inversion symmetry on it. Since [X y (0)] Ny = 
W y (0) [Eq. pO] )] and W y (0) is fixed by inversion symme- 
try to ±1, we have 



In Sec. IV B , we make the gauge choice with e l ® y ( x,ky " ) 
independent from X. We now drop the X argument and 
-1 and N y even; (94) write e i* v (k v ) m xhis enaD l e s us to further simplify the 
-1 and N y odd. e z ®v( k y) factors in Eq. (99). The starting point is the 

recursive definition of ^y^ k v) in Eq. (63). Define the 
shorthand notation 



[A„(0)] 



2iV„ 



1. 



(95) 



As shown in Appx. |C 2[ the curvature fluctuations are 
also subject to the inversion symmetry, 



U y (ky)U y (-ky - 1) = 1. 



(96) 



Recall from Eq. (|59) that (to y ) N y = n«" U v {k). We can 
recombine the U y factors in [Yl^ y U v (k)] 2 into pairs, and 
get 

H) 2Ny = I- (97) 



Combining this with Eq. (95), we have 



[A„(0)< 



l2JV„ 



1. 



(98) 



D. Inversion Symmetry of the Wannier States 



Ay(ky) — Ay(0, ky)Uy(ky). 



Then Eq. (63) could be rewritten as 



J& y (ky + l)-i&y(k y ) = A^(0) UJy ^ 
Ay(ky) 



(100) 



(101) 



thanks to the periodicity of e l ® y ( ky ^ in k y . We would 



with e l 



now like to relate, by recursion of Eq. (101), e l ® y ( ky ^ 
a i<$>y(-k y ) regu ^ wou ld be a product of 

~koi to ku 



\ y (0)ujy/Ay(K) along the path connecting —n, y uu 
Thanks to the periodicity of e t ® y ( ky " ) in k y: we can focus 
on the first Brillouin zone k y £ [0 . . N y ). Since — k y < k yi 
successive application of Eq. ( 101 ) leads to 



J®y(0,ky)-®y(0,-ky) 



[X y (0)ojy] 2ky / J{Ay(n)Ay(-n-l). (102) 



In Appx. |C 3[ we find the action of the inversion oper- 
ator V on the Wannier states \X,k v ): 



e *o,k ve -i*v(x,k v ) v \ Xt k ^ 



e -i* v (-jc,-fc„) \_ x ,-k y ) 
e -i*»(-^-i -fe») \-X - 1, -k y ) otherwise. 



(99) 
if W x (k y ) = 1, 



The product of A y factors in the rewritten form on the 
second line shows the inversion symmetry structure ex- 
plicitly: from the inversion transformation of A y (k x ) in 
Eq. (87) and that of U y (k y ) in Eq. (96), we have 



Ay(ky)Ay(-ky ~l) = ^ < *V ' *V + 1 . 



(103) 
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Plugging this into Eq. (102), we get the relative phase 



between the Wannier states which are inversion partners 



[A,(0Kr*. (104) 

Therefore, the inversion operation on the Wannier states 
[Eq. (99)] can be simplified to 



e^°[A y (0K]- 2fe ^|X,fc y ) 

= [\-X,-ky) XW X {ky)-- 

1 \—X — 1,—ky) otherwise. 



(105) 



This concludes our analysis of the action of inversion on 
the one-body states. 



Here the indices (s, f) are defined conjugate to (s,r), 
in Eqs. (77) and ( [78] ). Recall that the total momen- 
tum of the N e particles in state |\I/;s,r) is given by 
(K X ,K U ) = [K x + sN e ,CK y + C(r-5 y )N e ] mod (N x ,N y ) 
[Eq. ([69])] . The total momentum of its inversion part- 
ner given in the above equation is (—K x , —K y ) mod 



(N x ,N y ), as anticipated and shown in Appx. C5 

To summarize, our construction of lattice wave func- 
tions indeed preserves on the many-body level the inver- 
sion symmetry of the host Chern insulator. In contrast, 
the gauge and the original construction in Ref. [15] does 
not preserve inversion symmetry unless W y (0) — 1, as 
shown in Appx. [G] 



E. Inversion Symmetry of the Many-Body States 



VI. NUMERICAL TESTS 



When the one-body model is inversion symmetric, we 
require the FCI many-body wave functions on lattice to 
inherit the inversion symmetry of the FQH states on a 
torus. First we check whether the inversion transforma- 
tion is compatible with the Wannier-LLL mapping on the 
index level. Under inversion, the LLL state index trans- 
forms by j — >• —j mod = N x N yi while as shown in 
Eq. (105), the Wannier state index transforms by 



(X,ky) 



(-X,-ky) 
(-X-l,-ky) 



if W X {ky) = 1, 

otherwise. 



(106) 



In Appx. C4 1 we show that the ID index j x ^ k v defined 
by Eq. (43) transforms by 



— » 



-J 



X,ky 

■X,k v 



1 



if W^(0) = 
otherwise. 



-1 and N y odd; 



(107) 

Note that for a given system, all the j x > k v indices obey 
the same transformation rule. The transformation in the 
first case is not a good FQH symmetry on the single- 
particle level. However, on the many-body level, a si- 
multaneous change j to — j — 1 for all particles is a good 
symmetry, implemented by VT^ m . 

Finally, we are ready to investigate the inversion sym- 
metry of the many-body lattic e sta te g, r)i at defined 
in Eq. ([48]). Plugging Eqs. ( fl05h and ( [107] ) into the 
Wannier-LLL mapping in Eq. (48), we can get the in- 
verted amplitudes on the FQH side. Then, we can ap- 



ply Eq. (80) and relate it to another state in the g-fold 
multiplet and map it back to the lattice. This identifies 
the inversion partner of the many-body lattice state and 
obtains the many-body matrix elements of the lattice op- 
erator V in the g-dimensional subspace. The details are 
presented in Appx. |C 5| while the gist is the following. 
The q lattice states s,r)i at come in inversion pairs: 



oc 



P|*;s,r)iat 
f|*;5,f-l)i at ifW*(0) = 
^5 ^)iat otherwise. 



-1 and N y odd; 



We perform the Wannier construction in the five 
fermionic lattice models studied previously at filling 
v = 1/3: the checkerboard lattice model, 26 the Haldane 
model on the honeycomb lattice^ a two-orbital model 
that resembles half (spin-up) of the mercury-telluride 
two-dimensional topological insulator,^ the Kagome lat- 
tice model with spin-orbit coupling between nearest 
neighbors,^ and the spin-polarized ruby lattice modelP^ 
For each model, we use the parameter set specified pre- 
viously^ and focus on the lattice sizes and aspect ra- 
tios where the system has a threefold quasi-degenerate 
ground state in numerics well separated from the excited 
states. These are all inversion symmetric models, for 
which we list the values of W x (0), W y (0), and C for each 
model in Table |T] 

We obtain the "Laughlin" FCI states on the lattice 
using the formalism described in this paper, i.e. replacing 
\V) in Eqs. ([l5j) and ([70| by the Laughlin FQH states on 
the torus in the continuum. The threefold Laughlin FQH 
states are obtained by numerically diagonalizing the FQH 
model Hamiltonian with only the first pseudopotential, 
V m = 5 m ,i. Details of the continuum FQH Hamiltonian 
can be found in Appx. [A] As detailed in Sec. Ill C the 
continuum torus for the FQH calculations has the same 
aspect ratio as the lattice system, and the twist angle 
of the five models are given in the last column of Table [T] 
The three Laughlin states are exact replicas except for a 



TABLE I: Values of the Chern number C, the Wilson loops 
Wic(0) and W y (0), and the twist angle of the corresponding 
FQH torus, for the five lattice models at the studied param- 
eters. 



Model 


C 


W x (0) 


Wy(0) 


e 


Checkerboard 


-1 


+1 


-1 


tt/2 


Haldane 


-1 


+1 


-1 


tt/3 


Kagome 


+1 


+1 


+1 


tt/3 


Ruby 


-1 


-1 


+1 


tt/3 


Two-orbital 


+1 


-1 


-1 


tt/2 
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shift in the total momentum n y . As shown in Sec. |II A[ 



they can be recombined into |\P;s,r) and then mapped 
to the lattice using the Wannier basis. We will refer to 
the resulting lattice states as the Laughlin FCI states, as 
opposed to the Laughlin FQH states in the continuum. 

For all of the five models, we find that the Laughlin 
FCI states are indeed translationally invariant and inver- 
sion symmetric. In contrast, if we follow the prescription 
in Ref. |T5j the resulting many-body states break the in- 
version symmetry for the checkerboard lattice model, the 
two-orbital model, and the Haldane model. As explained 
in Appx. [G| the issue responsible for this is that the Wil- 
son loop W y (0) equals —1 in these models at the studied 
parameter set. 



A. Overlap 

The relevance of the Laughlin FCI states can be 
checked by computing the overlap with the ground states 
obtained from exact diagonalization of the FCI models at 
filling z/ = 1/3.™. 

The number and momenta of Laughlin FCI states 
(shown in Sec. |IVD ) match those of the exact ground 
states obtained by exact diagonalization. When both N x 
and N y are divisible by 3, all of the threefold states have 
total momentum zero. In such cases with more than 
one state in a certain momentum sector, the states can 
mix and we do not have a one-one correspondence be- 
tween the two groups of states to compute the overlaps. 
Instead, we first build the projector into the subspace 
spanned by the exact diagonalization states, then the 
projector into the subspace spanned by the Laughlin FCI 
states, and finally, we define the overlap to be the trace 
of the product of the two projectors, divided by the di- 
mension of the subspace, i.e. 



lat 1 



.<ED;i|*;s,r> 



lat 



(109) 



where |ED;z)i at , i — 1,2,3 are the threefold ground 
states obtained from exact diagonalization, and the in- 
dices (s, r) are summed over the threefold Laughlin FCI 
states in the same momentum sector. This naturally gen- 
eralizes the usual definition of overlap as the absolute 
square of the inner product. When the Laughlin FCI 
states appear at different momenta, we take the usual 
overlap with the exact diagonalization states. 

Shown in Fig. [6] against the Hilbert space dimension are 
the overlap values between the Laughlin FCI states ob- 
tained by the Wannier construction and the exact ground 
states at filling v = 1/3, of the five models on a lattice of 
size 4 x 3, 5 x 3, 6 x 3, 7 x 3, and 6x4. The only excep- 
tion is the two-orbital model: it does not have a gapped 
threefold quasi- degenerate ground state when N y = 3. 7 
At 6 x 3 the threefold states have the same total momen- 
tum, and thus there is only one overlap value for each 
model. At other lattice sizes, two of the threefold states 
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FIG. 6: Overlap between variational states and the exact 
ground states. The size of the lattice is labeled by N x x N y at 
top. The plot includes the checkerboard lattice model (pur- 
ple crosses), the Haldane model (red hexagons), the Kagome 
lattice model (green triangles), the ruby lattice model (blue 
Y-shapes), and the two-orbital model (cyan squares). The 
Wannier construction is performed using the \X,k y ) basis. 
At 6 x 3, the three states in a multiplet are in the same mo- 
mentum sector, and thus there is one overlap value for each 
model. At other lattice sizes, two out of the three states have 
identical overlap due to inversion symmetry. Also shown is 
the overlap between the FQH ground state of the Coulomb 
interaction at filling v — 1/3 and the Laughlin model state, 
denoted by black circles and labeled by the number of flux . 
As a baseline, we also plot the top 10 -6 and 10 -12 percentiles 
(gray lines) of the overlap values between random unit vectors 
as a function of the Hilbert space dimension. The overlap axis 
range is set to [0, 1] to contrast with Fig. [7] 



form an inversion pair and have the same overlap value. 
Also shown in Fig. [6] are the overlap values between the 
Laughlin state and the FQH ground state of the Coulomb 
interaction on a torus with unity aspect ratio, pierced by 
= 15, 18, 21, 24, 27 fluxes. 
We find that at comparable Hilbert space sizes, the 
overlap for the Kagome lattice model and the ruby lat- 
tice model are comparable with or even higher than the 
overlap for the Coulomb interaction in the continuum, 
while the overlap for the Haldane model, the two-orbital 
model, and the checkerboard lattice model have slightly 
lower overlap than the Coulomb interaction. We note 
one possible reason for the higher overlaps in the FCI 
systems. The lattice interactions used in the FCI models 
are Hubbard interactions between the nearest neighbors. 
Their short-ranged nature is preserved by the mapping 
to the LLL via the localized Wannier states, and thus 
the FCI interaction in the Wannier/LLL basis maybe 
be closer to the Laughlin pseudopotential Hamiltonian, 
compared with the long-range Coulomb interaction in the 
continuum. 

We emphasize that the overlap with the Laughlin FCI 
state should be interpreted in the context of the size of 
the Hilbert space. In Appx. [H] we find that the proba- 
bility for two random complex unit vectors in a space of 
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FIG. 7: Overlaps between variational states and the exact 
ground states using the prescription in Ref. 1151 Inversion 
breaking for the checkerboard lattice model and the Haldane 
model is clearly visible at 4 x 3, and present (although not dis- 
cernible in this figure) at other lattice sizes. Inversion break- 
ing for these models is also present at 7 x 3 and 6x4, even 
though obscured by the finite resolution of the plot. The in- 
version symmetry of the two-orbital model at 6 x 4 is also 
broken. 



dimension d to have an overlap larger than a given value 
is exponentially small at large d. This is illustrated in 
Fig. [6] by the plot of the top 10 -6 and 10 -12 percentile 
values of the overlap between random vectors, as a func- 
tion of the Hilbert space dimension. 

We observe a slight decreasing trend in the overlap as 
the system size increases, much slower than the exponen- 
tial decay of the random overlap. There is also an upward 
kink when going from 7 x 3 to 6 x 4 for the checkerboard 
lattice and the Haldane models. This agrees with the 
finding in Ref. 4 that the topological phase is more sta- 
ble when the aspect ratio is closer to unity. 

We have also tested the performance of the proposal 
in Ref. [15] While the original proposal was given in the 
continuum limit N Xl N y — >> oo (formulated using deriva- 
tives such as (k xi ky\dk x \k xi k y )) and on a cylinder, as 
detailed in Appx. [G| we extend this proposal to handle 
finite-size lattice models with multiple sublattices and 
periodic boundaries, but do not change the gauge used 
there. We repeat the overlap calculations using this al- 
ternative prescription. The results are compiled in Fig.[7[ 
The overlap values for the Kagome lattice model and the 
ruby lattice model are comparable with those calculated 
using our proposal. However, the overlap values for the 
checkerboard lattice model, the two-orbital model (which 
is the model in the original discussion 15 ), and the Hal- 
dane model are much lower than the values obtained from 
our prescription. At 6 x 4, the overlap values for these 
three models fall below 0.04. Also, the constructed many- 
body states break the inversion symmetry in these mod- 
els. The culprit is explained in Appx. |F3| and Appx. [G] 
In short, when W y (0) ^ 1, the gauge choice a y = used 
in Ref. [15] does not mimic the behavior of the LLL or- 




FIG. 8: Overlaps between the Laughlin FCI states and the 
exact diagonalization states, as a function of the shift pa- 
rameter S y in the Wannier-LLL mapping. Please refer to the 
legend and caption of Fig.[6]for the annotation of scatter sym- 
bols. The left (right) panel shows the calculations on a 6 x 4 
(6 x 3) lattice. For clarity, for each model at each S y we show 
the average overlap of the threefold states. The two-orbital 
model is not included in the 6x3 calculation, as it does not 
have a gapped topological ground state at this lattice size. 
The values of 5 y that produce significant overlaps are in full 
agreement with Eq. (94). 



bitals: enforcing a y = in the interior of the Brillouin 
zone mandates that the connection on the Brillouin zone 
boundary has to take care of the gauge- invariant Wil- 
son loop, and deviates from a y = 0. The prescription 
in Ref. [15] does not properly handle the cases where the 
Wilson loop W y (0) is not unity, which include the above 
three lattice models at the studied parameters, as shown 
in Table [I] 



B. Mapping Parameters Revisited 

In the previous section, the shift parameter S y for the 
Wannier-LLL mapping and the twist angle of the con- 
tinuum torus used in the FQH calculations were set ac- 
cording to the discussions in Sec. |III C| We now revisit 
the choice of S y and by checking explicitly how it affects 
the overlaps with the FCI ground states. 

First we perform the Wannier constructions for the five 
models using the angle specified in Table [IJ but we vary 
the shift parameter S y over [0 .. N y ). In Fig. [8j we show 
the dependence of the overlaps on S y for 6 x 4 and 6x3 
lattices. At N y = 4 (even), for the three models with 
W^c(0) = +1? namely (see Table [I]), the Kagome lattice 
model, the checkerboard lattice model, and the Haldane 
model, the Laughlin FCI states have significant overlaps 
with the exact diagonalization ground states only at S y = 
0; for the two models with W x (0) = —1, namely, the ruby 
lattice model and the two-orbital model, the significant 
overlaps appear at S y = 2 = N y /2. At N y = 3 (odd), the 
significant overlaps for the three models with W^O) = 



+1 are still at 5, 



0, 



but for the ruby lattice models 
with Wa;(0) = — 1, the peak in overlap is shifted to S y = 
1 = (N y — l)/2. 65 The above results are in full agreement 
with the choice of S y in Eq. (94). 



We now put S y to the values specified by Eq. ( 94 ) , but 



18 



1.0 
0.8 



& 0.6 

> 

O 0.4 



0.2 
0.0 















Haldane 

Checkerboard 

Ruby 

Kagome 

Two-orbital 







tt/3 



tt/2 




FIG. 9: Overlaps between the Laughlin FCI states and the 
exact diagonalization states on a 6 x 4 lattice, as a function 
of the twist angle of the continuum torus for the FQH cal- 
culations. For clarity, for each model at each 5 y we show the 
average overlap of the threefold states. For each model, a 
peak in overlap is clearly visible at the value of 6 given in 
Table [I] in accordance with the discussion in Sec. |III C] 
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FIG. 10: Overlap optimization for the Laughlin lattice states 
by tuning the single-particle gauge. The calculations are per- 
formed on a 6 x 4 lattice. Please refer to the legend and cap- 
tion of Fig.[6]for the annotatio n of scatter symbols. Shown in 
a) is the overlap upper bound ( 110 ) obtained by removing the 



relative phase at each component between the exact ground 
state and the Laughlin lattice state, against the actual over- 
lap. Shown in b) is the supremum (the least upper bound) 
of the overlap, obtained by a brute force optimization of the 
single-particle gauge of the Bloch states, against the actual 
overlap. 



vary the twisted angle of the continuum torus for FQH 
calculations. The overlaps as a function of are shown 
in Fig. [9] For each of the five models, we find a clear 
peak in overlap centered around the value of suggested 
in Sec. |III C\ as shown in Table [TJ This confirms that our 
choice of the twist angle is appropriate. 



C. Gauge Optimization 



As detailed in Sec. IV B the phase choice e 



of the Wannier basis is based on the similarity with the 
LLL orbital behavior. There are other possible alterna- 



tive gauge choices. We now try to understand how close 
our gauge fixing procedure gets to several brute-force op- 
timized gauges. 

We first examine the limit where the phases of the 
Slater determinant weights in the variational trial states 
are tuned to match perfectly with the phases in the ex- 
act diagonalization ground states. Eq. ( pfp] ) shows that 
the absolute value of each component ({k x , k y }\^; s, r)\ at 
of the many-body state in the Bloch basis does not de- 
pend on the parallel-transport gauge in the k x direction, 
nor on the gauge choice of the Wannier states if we limit 
ourselves to a set of e l ^ v ^ x ' ky ^ phases that have no real 
X dependence. Inspired by this, we remove the relative 
phase from the Bloch-basis overlap formula between the 
Laughlin lattice state |^;s,r)i at and the corresponding 
exact diagonalization ground state |ED;s,r), and calcu- 
late 



£ . 

{k x iky } 



iat(ED; s,r\{k x ,k y })({k x ,ky}\^; s,r) ]at 



■ (no) 



This real number serves as an upper bound to the ac- 
tual overlap. It neglects the phases of the Slater deter- 
minant weights in the many-body wave functions.^ In 
Fig. 10 1), the numerical values of this upper bound are 
shown against the overlap values from our prescription 
for the five lattice models. We find that the upper bound 
is still far from unity for the three under-performing mod- 
els, namely, the checkerboard lattice model, the Hal- 
dane model, and the two-orbital model. This can be 
interpreted in two possible ways: The X dependence of 
e i<5>y(x,k y ) ma y -fo e essen tial for the three models, or, the 
ground states of the three models may be simply more 
complicated than the Laughlin lattice states. Our results 
are not sufficient to distinguish between the two scenar- 
ios. 

The upper bound is higher than the actual overlaps 
from our prescription, which would naively suggest that 
our prescription is not optimal. Thi s is n ot true. To reach 
the upper bound calculated in Eq. ( 110), we need to tune 



the phases of all the Slater determinant weights in the 
many-body state. The number of such phases is on the 
order of (N x N y )- choose- N e . We in no way have the free- 
dom to tune these phases independently. The most we 
are allowed to do is tune the phases of the N x N y Wannier 
states independently. The least upper bound, i.e. the 
supremum, is hence found by tuning the single-particle 
gauge to reach maximum overlap. In general, this vari- 
atio nal o verlap value cannot reach the upper bound in 
Eq. (fTTob- 



The actual supremum can be found by a brute-force 
optimization. We introduce N x x N y variational gauge 
degrees of freedom e lTkx > k y for the Bloch states \k x ,k y ). 
This transforms the Bloch amplitudes by 



({fca.,fc y }|*;s,r)iat 



' ky ({k x ,k y }\^;s,r)i at . 



We optimize the average overlap with the threefold ex- 
act ground states by tuning these N x N y parameters. As 
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explained in Sec. IIIB 



3 zTo > k y covers the degrees of free- 
dom carried by the Wannier phase choice e l ® y ( ky \ We are 
not able to perform the optimization on e L ^ y ^ x,ky ^ with 
non-trivial X-dependence due to the prohibiting compu- 
tational complexity. We use the Nelder-Mead (downhill 
simplex) algorithm 3 -^ i Q fi nc [ the local overlap maximum 
around our prescription. The convergence threshold is 
10 -4 on the overlap (output) and the phase angles of the 
gauge transform (input). 

In Fig. [10)3), the suprema are shown against the 
overlap values from our prescription. Comparing with 
Fig.[lQ^), we find that the supremum is much lower than 
the upper bound (110) for the three under-performing 



models. This reduction is not surprising, since the varia- 
tional degrees of freedom are far fewer than the number 
of many-body amplitudes. The proximity of the data 
points to the diagonal line in Fig. 10 d) suggests that our 



prescription is very close to the optimal one if we refrain 
from introducing X-dependent phases in \X,k y ). From 



the values of e 



at the supremum, we find that to 



improve our prescription, we need to move slightly away 
from the parallel transport gauge [Sec. IIIB in the k x di- 
rection. This may be attributed to the curvature fluctu- 
ations. We leave explorations in this direction for future 
work. 



D. Entanglement Spectrum 

Overla p only provides an overall, sometimes inaccu- 
rate^^ estimate on the similarity between two states. 
To compare the correlations in the Laughlin FCI states 
and the ground states, we now turn to their particle en- 
tanglement spectra. 

First introduced by Li and Haldanef^ the entangle- 
ment spectrum examines the information imprinted in 
a many-body state. By effectively performing a spectral 
decomposition of the many-body state, the entanglement 
spectrum reveals the excitations supported by the state. 
Of particular interest in our study is the particle entan- 
glement spectrum (PES), which encodes the characteris- 
tics of the quasihole excitations from the ground state. 33 
We divide the N e particles into two groups A and B of 
Na and Nb particles respectively, trace out the degrees 
of freedom carried by the particles in £?, and examine of 
the negative logarithm of the eigenvalues of the reduced 
density matrix qa> 

The Laughlin FQH state exhibits a highly non-trivial 
pattern in the PESP^MJi n each momentum sector, the 
number of entanglement levels at finite entanglement en- 
ergy matches the counting of the Laughlin quasiholes of 
Na particles in qN e fluxes. The total number of levels, 
i.e. the number of non-zero modes of the reduced den- 
sity matrix, is much smaller than the dimension of the 
reduced density matrix. 

The PES of various FCI ground states at filling v = 1/3 
has a low entanglement energy structure similar to the 
PES of the Laughlin state PEGS it also displays a non- 
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FIG. 11: Particle entanglement spectrum of the Kagome 
lattice model of TV = 8 particles on the N x x N y = 6x4 lattice. 
Shown in a) and b) are the spectra of the Laughlin FCI state 
and the exact ground state, respectively, with Na = 3. Shown 
in c) and d) are the corresponding spectra with Na = 4. In 
each momentum sector, we calculate the number of admissible 
configurations n, and color the lowest n entanglement energy 
levels in red and the rest in blue. The structure shown in the 
left panel exhibits the defining characteristic of the Laughlin 
state. 



universal high entanglement energy structure separated 
by an entanglement gap from the Laughlin-like struc- 
ture. The number of levels below the gap at each total 
momentum matches the number of quasiholes folded to 
the lattice Brillouin zoneP An ideal Laughlin state in the 
FCI would exhibit an infinite entanglement gap. 

For the Laughlin FCI model state obtained from the 
Wannier construction, we find a large number of zero 
modes in the reduced density matrix - corresponding to 
an infinite entanglement gap. The number of the finite 
entanglement energy levels is given exactly by the count- 
ing of Laughlin quasiholes folded to the lattice Brillouin 
zone. Writing the Laughlin FQH state in the Wannier 
basis is just a unitary transformation. Thus it does not 
modify the rank of the reduced density matrix, and it pre- 
serves the entanglement spectrum of the Laughlin FQH 
state. An example is shown in Fig. 11 The entangle- 



ment energy levels below the gap of the exact diagonal- 
ization ground state exhibit almost identical structure in 
the spectra of the Laughlin FCI state and the ground 
state. The Laughlin FCI state captures a very important 
feature of the correlations in the FCI ground states. This 
further corroborates the proposal^ that the quasihole 
counting in the PES can be interpreted as a signature of 
a FQH-type of topological phase. 



VII. ISOTROPY/NEMATICITY OF THE 
WANNIER CONSTRUCTION 



An undesirable feature of the Wannier construction is 
that it seems to explicitly break the isotropy between the 
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FIG. 12: Single-particle orbitals in the lowest Landau level 
(LLL) and the Chern band in the alternative gauge choice. 
In the left panel we show the Landau orbitals in the LLL on 
a torus. The two fundamental cycles are marked by Li and 
L12, and the twist angle is labeled by 0. The Landau orbitals 
are plane waves in the e v direction that are localized in the 
direction perpendicular to propagation. In the right panel we 
show a ID Wannier state localized in the b2 direction and 
in the lowest band of the Kagome lattice model. The size of 
the spheres depicts the weights of the Wannier state on each 
lattice site. 

two axes of the lattice system when the original Hamil- 
tonian is not fourfold rotationally symmetric.^ The for- 
malism singles out one direction in the 2D lattice along 
which the single-particle states (both on lattice and in 
the LLL) localize. The two alternative choices of localiz- 
ing along the x/y direction are not guaranteed to produce 
the same set of g-fold many-body states for the FCI. 

In the following, we show that the weights of the cor- 
responding Bloch-basis Slater determinants in the two 
sets of many-body states constructed from the two al- 
ternatives have the same absolute value, but the phases 
may differ. However, in the limit of flat curvature, the 
two alternatives produce exactly the same set of g-fold 
many-body states. (In this paper we do not explore other 
possible localization directions of the Wannier basis.) 

So far we have been using single-particle states \X,k y ) 

I 



The state \j)' v is localized along the line y = xcot0 + 
jL 2 /N ( f > . The center line moves in the positive e y direc- 
tion when we increase j. 



The Slater determinant basis for the 7V e -electron 

Hilbert space generates a representation of the many- 



that are plane waves in the e y direction. We now exam- 
ine the alternative choice of using single-particle states 
\Y,k x ) that are localized in the e y direction, and the cor- 
responding alternative basis in the LLL. 



A. FQH States in a Rotated Gauge 

We look for a one-body basis in the LLL that are lo- 
calized in the e y direction, as illustrated in Fig. [12] The 
"natural" gauge in this case is aligned along the e v direc- 
tion. We setup a rotated coordinate system (u,v) with 
the v axis pointing in the direction of e v : 

u(x, y) = x cos — y sin 0, 

„ L (in) 

v(x, y) = x sin + y cos 0, 
where is the angle between the two fundamental cycles 
of the torus. Then the new gauge is the Landau gauge in 
the rotated coordinates A!(x,y) = Bu(x,y)e v . This is 
related to the original Landau gauge A = Bxe y by the 
gauge transform A' — A = V%, with V = (c%, d y ) and 

X&y) = — ^{x,y)v{x,y) - xy\. (112) 

Since the guiding-center momentum is gauge- 
covariant, the guiding-center translation operator in the 
primed gauge is given by 

Ifta) = e iex ' h Ti(aL)e- iex / h . (113) 

From this we can obtain the primed variants of the many- 
body translation operators discussed earlier. 

We denote by \j)' v and the single-particle states in 
this gauge that satisfy the guiding-center periodic bound- 
ary conditions 7^(Li) = T((L 2 ) = 1. The prime high- 
lights the alternative gauge choice A 7 , while the super- 
script v emphasizes that this state is propagating in the 
e v direction (and localized in the e y direction). The wave 
functions of these states resemble those in the unprimed 
gauge [Eq. (g], 



e -lu(x 1 y)] 2 /(^%) m (H4) 
I 

body translation operators. In particular, we find 

^i{i}); = e - i27rEj/N i{i});, 
T%MY v = \{j -?»;> (115) 

T^\{j})' v = \{j + l})' v . 



n u(x, y) + iv(x, y) tt Loe %e . , T n9 
-27r(j + nN^) V ,y> V ,y> + + nN^) 2 
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In contrast to the unprimed transformations in Eq. ([9|, 
the translations in the Li (rather than L2) direction are 
diagonal, and the eigenvalues of the diagonal operators 
differ by a complex conjugate. 

We can diagonalize the interaction system with the 
help of the many-body translational symmetries. For ob- 
vious reasons the eigenstates are just the gauge-transform 
of the states discussed in Sec. II A Recall the state 
diagonal in (T^Tj^) with wave number (n xi n y ). De- 
note by \^f)' the gauge-transformed state. Then we have 



i2irn x /N 



= e - i27iK y' N *\^)' . 



(116) 



Slater determinant state formed by a Fourier series of 
the Landau orbitals \m) transforms in the same way as 
under the many-body translation operators. We 
can actually prove that 



\j)v 



i(f-f) 



-i2irjm/N ( ( ) 



m). 



(120) 



This behavior under a Fourier transform is expected of 
the shifted Gaussian form of the single-particle wave 
function. The details are presented in Appx. [D] 
B. Alternative Wannier Construction 



And the simultaneous eigenstates of S' x , R' y , and the 
Hamiltonian can be constructed from \^)' in a parallel 



manner to the construction of s, r) state in Sec. II A 



The result is of course given by the gauge transform 
|*;*,r)'. 

We can express the \Sfr;s,r}' states by the amplitudes 
s, r)' . We emphasize that 



l{{jms,r)'^({j}\y;s,r), 



(117) 



because the single-particle state \j)' v is not the gauge 
transform of the state \j). We now seek a link between 
the two unequal amplitudes in Eq. ( |117[ ). First, we can 
transform | j)' v back to the unprimed gauge, using y) 
given by Eq. \112\ . The resulting states \j) v are given 
by the wave function 

<t>)(x,y) = e^^-^'^^'^/^ 2 ^)^^,^). (118) 

These state are still propagating in the e v direction, as 
signified by the super /sub-script v. We thus arrive at 
an alternative single-particle basis \j) v of the LLL in the 
Landau gauge. 

Thanks to the invariance of the inner product under 
the gauge transform, the many-body amplitudes satisfy 



l({jms,ry = v {{j}\*;s,r). 



(119) 



The last ingredient is the expansion of the state \j) v in 
the \j) basis. We have already seen some hint from the 
two representations of the many-b ody translation opera- 
tors. Comparing Eq. (pi and Eq. (115), we find that the 



The Wannier basis |Y, k x ) 



localized in the e y direction 



(see Fig. 12) can be obtained trivially from |X, k y ) by 



exchanging all x and y and replacing C by — C. Using this 
substitution rule, we can define 5 X1 uo x1 U x (k x ) alongside 
their sub-?/ counterparts. 

We link the \Y,k x ) state to the LLL orbital \j Y > kx )' v . 
The index mapping is given by 



Jv 



YN X - Ck x + 5 X 



(121) 



Here we put a subscript v to distinguish it from the j x > k v 
mapping defined earlier. Then we can transcribe the 
FQH states \^;s,r) f to the lattice. The resulting many- 
body lattice states |\I/;s,r)( at are translationally invari- 
ant. Despite the name, the new lattice states |\£;s,r)( at 
are not the gauge transform of the states s, r)i at , since 
they are constructed from different Wannier mappings 
not related by any gauge transform on a lattice at zero 
magnetic field. We emphasize that only in the LLL does 
the prime denote the rotated Landau gauge. 

The total momenta of s, r)[ at are given by 



K x = C[k x + (s + S x )N e ] mod N x , 



K y + rN e mod N y . 



(122) 



The counting of the q states in each total momentum 
sector is exactly the same as the counting of the unprimed 
lattice states. The amplitudes in the Bloch basis are 
given by 



({k x ,k y ms : r)[ at = l[- 

By comparison, recall that 
({k xi k y }\^;s,r)\ ax = 



ui x u x (K)u: x A x ^o)u: y A y (k x ,K)) jw % 



(123) 



y {y} 



[Ay(0)] fc " [A,(fc,)] fc - 
Yt v U y {n) U k K v Ay(0, k) Yt A x (k, k y ) 



{X} 



J70 



22 



C. Lattice Amplitudes 

Now we are ready to check the isotropy/nematicity of 
the lattice states obtained by the Wannier constructions. 

To avoid the unnecessary clutter, here we focus on the 
simplest case with C = +1 and 5 X = S y = 0. The generic 
case is treated in Appx. [Eland the conclusion is the same. 

In this case, the Wannier mappings are reduced to 



jX'kv =XNy+ky, 



Jv 



YN x -k x . (124) 



Recall from Sec. |IIIC| the subtle point that the above 
mapping holds only if (k x ,k y ) is in the principal Brillouin 
zone, in this case (—N x .. 0] x [0 .. N y ). To minimize the 
trouble, we set all the N e mome nta in {h x , k y } on the left 
hand side of Eqs. ( pfo] ) and |123 to the principal Brillouin 



zone. 

Note that the fractional factors enclosed by the curly 
braces in those equations are invariant under a shift of 
Brillouin zones: for example, (uj x ) kx / Yl^ x U x (k) does 
not change if we shift k x — >• k x + N x . Therefore, we 
can always pull (k x ,k y ) back to the first Brillouin zone 
[0 .. N x ) x [0 .. N y ) when evaluating the prefactors in the 
curly braces. 

Thanks to C = +1, S x — S y = 0, the dependence 
of (K X: K y ) on (s,r) is identical for both s,r)i a t and 
s, r)[ at . Hence we try to relate the s 7 , r')' and 

(0"}|*;s,r) with 



( 5 7 ,r 7 ) = ( S ,r). 



(125) 



We first plug the Fourier transform in Eq. (120) into Eq. (119). We can rewrite each sum over m in the Fourier 
transform as a double sum over all unit cells X and all l y points in the principal Brillouin zone [0, N y ). We find that 



^2 e -i2nj:k x X/N x ^ e i27v ^2 l y Y/N ye -i27r ^2 k^y/N^ ({j X ^ }\^f • S , r) , (126) 



-N € 

</> {X} 



{ly} 



where e %1 = e zNe d 2) is a constant phase factor. Here we have used e l2 ^ XYN xN y /N (j) _ thanks to N x N y = 
for \C\ = 1. Therefore, we have 



^1 



1 ^2 e~ i2n ^ k v Y / N v v ({ j Y > kx }|#; s 7 , r')' = e i7 e~ i27T ^ k * k y/ N <t> 1 ^ £ k x x/N x ({j*,k y } |^ ; 5? 

s v\ 



y {Y} 



x {X} 



(127) 

Here the sum over each Y produces a v . Since both k y and l y are limited to [0 .. N y ), this is just Sk y ,i y and 

thus we can remove the sum over {l y }. Putting the above equation into Eq. (123), we find that 



{{k x ,k y }\^s',r')[ &t = e^\{ 



K) fc * [A,(0)] fc * l\(k x )P e -i2,ck x k v /N, ' 

It u x (k) l£ x M*,o) Yt y A y (k x , K ) 



_J_^ e -i2^Y.k,x/N x ({j-«.*»}|* ;S)r ). (128) 
/N' e s v\ 



x {X} 



The above equation holds (up to an inconsequential 
overall constant phase factor) for generic (C,5 x ,5 y ) val- 
ues, as proved in Appx. [EJ 

Disregarding the overall constant phase e* 7 , the 
only difference between ({k x , k y }\^; s 7 , r f )[ at and 
({k x , k y }\ty; s,r) \ a t [Eq. ( |70| )] is the phase enclosed by 
the curly braces. Therefore, the two alternative choices 
of the Wannier localization direction produce two sets of 
many-body lattice states with exactly the same absolute 
value for each Bloch-basis component, 



({k x ,k y }\V;s f ,r f )[ at = <{fcs,fcj|tf;s,r)i at . (129) 



notations for the phase prefactors in Eqs. (fTOI) and (123): 



^ = K)^ [A,(0)]^ [X X (ky 



^4>K x ,K y 



g K X , K y 



Yt v U y {n) Yt v A y (0, K ) Yt x A x (K,k v ) 

{u x ) k * [\ x (0)] k * [X y (k x )P 
It U x (k) \t M", 0) U k K y Mkz, «) 



(130) 

Intuitively, e l(pkx > k y comes from the parallel transport in 
the Brillouin zone along 

(0,0) <- (0,^) <r- (k x ,ky), 

while the phase e l ^ kx ' k y comes from the parallel transport 
along 



We now examine the phase part. Define the shorthand 



(0,0) <- (k x ,0) <r- (k x ,k y ). 
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FIG. 13: Isotropy of the Wannier construction, measured by 
the overlap between the two sets of Laughlin lattice states con- 
structing using the \X,k y ) and \Y,k x ) Wannier bases. Please 
refer to the caption and legend of Fig. [5] for the annotation of 
scatter groups and lines. Inversion symmetry symmetry can 
be observed similar to Fig. [6] 



that the phase difference between ({k x , k y }\^; s, r)[ at and 
({k x , k y }\^f; s, r)i a t does not vanish in the continuum 
limit if the curvature is not flat. It cannot be eliminated 
by any point group symmetry, either. 

A detailed analysis of the impact on the isotropy of the 
Wannier construction is hard to carry out. But we can 
still examine the isotropy or nematicity numerically. We 
find strong isotropy in the five FCI models despite fluctu- 
ating Berry curvature. Shown in Fig. [13] is the isotropy 
measured by the overlap between the Laughlin lattice 
states constructed using the \X,k y ) and the |Y, k x ) Wan- 
nier bases. The isotropy is higher than 0.997 for the 
ruby lattice model and the two-orbital model at all sys- 
tem sizes checked. A decreasing trend is observable for 
the checkerboard lattice model and the Haldane model, 
but the slope is much slower than that of the overlap 
values with the exact ground states. 



VIII. DISCUSSION AND CONCLUSION 



As these phases are not identical for all Slater deter- 
minants, the Wannier wave functictions exhibit a certain 
degree of nematicity. However, in the limit of flat Berry 
curvature, we expect that the two different routes of par- 
allel transport picks up a relative phase of e %27vk x k y/ N ^ 
for C = 1, which is canceled exactly by the e ~' l2 ^ k xky/ N 4> 
phase from the basis transform from \ j) v to \j). 

The algebra indeed works out. As explained earlier, we 
can shift (k xi k y ) to [0 .. N x ) x [0 .. N y ) when evaluating 
the phase prefactor. For flat curvature, we have uj XjV = 



x,y 



1, and the Wilson loop around the rectangle is 



t 



[K, Ky) 

_ e i27rCk x k y /N (j) ^ (131) 

Recall that X x (k y ) and X y (k x ) are related to the Wilson 
loops W x (k y ) and W y (k x ) by Eq. (30), and thus 



X y (k x ) 


ky 


\ A*(0) "I 


[ A„(0) J 




_X X (ky) _ 



k: x 



p i2-KCk x k y I 



Putting all these together, we find as claimed 

i27rCk x ky/N (f) 



Ky 6 



(132) 



(133) 



and thus the Wannier construction is indeed isotropic in 
the limit of flat Berry curvature: 



laf 



(134) 



D. Numerical Check 



Realistic lattice models inevitably introduce fluctu- 
ations in the Berry curvature. The fluctuations per- 
sist in the thermodynamic limit N e — >• oo. We note 



In this paper, we have studied the Wannier construc- 
tion of FQH states for fractional Chern insulators. De- 
spite being first proposed more than a year agcP^, no 
overlap studies with the exact diagonalization have been 
reported. If one did the comparison naively, low over- 
lap would be obtained (see Fig. [7]). We have proposed a 
new, related prescription which fixes several outstanding 
issues in Ref. 15, In particular, we highlight the gauge 
freedom in the ID localized Wannier states and the re- 
striction imposed by the Wilson loops. The key point is 
that the relative phase between adjacent Wannier states 
must follow the same pattern as the LLL orbitals in the 
Landau gauge, otherwise the states obtained from the 
FQH-FCI mapping are not relevant to the FCI problem. 

We describe and justify in details a finite-size prescrip- 
tion for a generic lattice model with multiple sublattices. 
We provide an explicit, step-by-step recipe to construct 
the counterpart of any FQH wave function on a Chern 
insulator. Our prescription is tailored for the torus geom- 
etry to make contact with existing numerical studies. It 
preserves the full translational invariance in both direc- 
tions, as well as the inversion symmetry. Moreover, in the 
limit of flat Berry curvature, the constructed many-body 
states are the same using the Wannier basis localized in 
either x or y direction. 

We find the major obstacle to the Wannier construc- 
tion to be the fluctuations in the Berry curvature. We try 
to accommodate the fluctuations in the localization direc- 
tion in the gauge choice of the Wannier basis. The fluc- 
tuations in the other direction keeps the Wannier centers 
from evenly spaced like the LLL orbitals, and we find no 
way to ameliorates the situation without spoiling trans- 
lational invariance. We have no real control over the ne- 
maticity developed from the curvature fluctuations. This 
complication from curvature is an important feature that 
sets the FCI apart from the FQH in the continuum. 

In an earlier paper^, two of us discovered a folding rule 
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to relate the counting of the FQH states on a torus to the 
counting of the FCI states at each total momentum. It 
was established by arguing that the FCI phase on a lat- 
tice is in the same universality class as the FQH phase 
on a torus with ^-function lattice pinning potentials. In 
this paper, the Wannier construction reorganizes the de- 
generate FQH model states and produces a set of model 
lattice states. These states have exactly the same count- 
ing in each total momentum sector of the energy and the 
entanglement spectra as predicted by the counting rule 
in Ref. 9 , This provides a concrete implementation of the 
folding rule that connects the FCI and the FQH effects. 

We perform the Wannier construction for the Laugh- 
lin state numerically in the five lattice models known to 
support a fractionalized phase. We find that the model 
state obtained from the Wannier construction and the 
actual FCI ground state have consistently large overlap 
in all cases, and their entanglement spectra exhibit very 
similar structures. The many-body physics in Chern in- 
sulators at filling v = 1/3 is indeed a close parallel to the 
Laughlin-type of FQH physics in the continuum with a 
strong magnetic field. 

Using these five examples, we demonstrate that our 
prescription for the single-particle gauge is at or very 
close to the local optimum in terms of the overlap with 
the exact ground states. We also show that the our pre- 
scription has strong isotropy by checking the overlap be- 
tween the many-body states constructed using the Wan- 
nier bases localized in either x or y direction. 

Our results provide another comparison across the ar- 
ray of known FCI models. Among the five models that 
we have checked, the ground states of the Kagome lattice 
model and the ruby lattice model have the largest overlap 
with the Laughlin model state. This is consistent with 
the finding in Ref. 7 that these two models host the most 
stable Laughlin-type FCI phase among the five. 

We note a few interesting directions for future work. 
Two of us recently studied the thin-torus limit of the 
FClpy j n ii m it, the projected interaction is domi- 
nated by classical electrostatic terms and allows an exact 
or perturbative solution. This may provide a better un- 
derstanding of the varied performance of the Wannier 
construction in different models. 

In this paper we focus the numerical tests on the 
Laughlin state. The general formalism of the Wannier 
construction, however, is able to handles more compli- 
cated FQH states. Lattice analogs of composite fermion 
states^ and non-Abelian FQH states such as Moore- 
ReacP^ and Read-Rezayi stated^ have been identified 
in a few Chern insulators! 7 * 9 * 23 * 3 ^ Moreover, it was pro- 
posed that topological flat bands with Chern number 
|C| > 2p3EI ma y be may be understood as a multi- 
component FQH systemP^ It would be interesting to see 
how well these states could be captured by some estab- 
lished FQH trial wave functions mapped to the lattice by 
the Wannier construction. 

Time-reversal invariant fractional Z2 topological insu- 
lators have been constructed from two copies of FCI and 



the field theory descr ibing their low-energy excitations 
have been established ! 44 * 45 * More recently, the algebraic 
structure of the projected density operators for topologi- 
cal insulators in three^or higher dimensions 4 ^ have been 
studied. With necessary generalizations, the gauge-fixed 
Wannier construction may be useful for the study of the 
wave functions of these more complicated systems. In 
particular, the case of 2D FTI (quantum spin Hall) can 
be tackled similar to the FCI problem, as the FTI system 
becomes decoupled in the Wannier basis. 

Another interesting approach to build tr ial wa ve func- 
tions for the FCI is the parton constructio n 48 * 49 !, initially 
pioneered by We n 150 * 51 ! to construct the FQH wave func- 
tions in the continuum. This alternative is plagued by a 
critical issue at the moment. To achieve fractional filling 
v = 1/q of electrons, the number of partons in a fully 
filled band should be reduced accordingly. For the orig- 
inal FQH system, this is not a problem since the size 
of the parton LLL orbital increases automatically due 
to the decrease of charge. On a lattice, this change has 
to be put in by hand through a unit-cell enlargement. 
Refs. 48 49] discussed how to preserve the lattice symme- 
tries in the electronic wave functions, but their solutions 
involve partonic bands that are drastically different from 
the original electronic band: even the band topologies 
may be totally different. It would be interesting to find 
a cure to this unsettling problem, to develop a concrete 
prescription for correct torus degeneracies, and to com- 
pare with the Wannier construction and the exact ground 
states. 
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Appendix A: FQH Projected Hamiltonian 

In this Appendix we give the expression of the FQH 
interacting Hamiltonian projected to the lowest Landau 
level (LLL) on a twisted torus. This is a trivial general- 
ization of the formula given in Ref. |52j 

In the limit of large cyclotron energy, we keep in the 
Hilbert space only the states in the LLL. Since the kinetic 
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energy is fully quenched, the Hamiltonian is simply the 
projected density-density interaction term V. 

Define the field operator ipj to annihilate the Landau 
orbital state \ j). The LLL-projected density operator can 
be written in momentum space as 

p q = Eo'ii e " iq ' ? i^)V'] 1 V'i 2 (Al) 

31,32 

d?e-^' ? <t>l (?)<M?)^i 2 , (A2) 

31 J 2 

where the reciprocal lattice vector q is given by q = 
(ZiGi + q2^2 with (#1,92) £ the continuous coordi- 
nate r is integrated over the principal region of the torus, 
and the orbital index ji, 22 are summed over [0 .. N<f>). 
Following the detailed steps laid out in Ref. [53j we can 
simplify the above formula and get 



"|q| 2 ^/4 



-i2Trq 1 (j+q 2 /2)/N (f ,-/A 



j Vj+q2 • 



(A3) 



In the second-quantized form, the projected two-body 
interaction can be written as 

^=^^£W q , (A4) 

where the reciprocal lattice vector q = q\G\ + #2^2 is 
summed over the infinite grid (^1,^2) £ Plugging in 
p q , we find 



V 



1 



2L1L2 sin# 

q ji j2 

p i27rq 1 (j 1 -j 2 -q2)/N (t> i\ j\ 



(A5) 



where Vq = J dre~* q ' r V(r) is the Fourier transform of 
the interaction potential V(r). 16 The Fourier coefficient 
can be expressed in terms of the Haldane pseudopo- 



tentials V m ^ 



F q = 4^|^K n L m (|q| 2 /|), 



(A6) 



where L m (x) are the Laguerre polynomials. 

As a side note, we point out an interes ting link betw een 
adjacent \j) orbit als. Combining Eqs. (Al) and (A3) at 
q = G2, we find 

{j\e- iG2 '*\j + 1) = e-l G2 l 2 ^/ 4 e R+. (A7) 

This can be viewed as the defining property of the par- 
ticular gauge choice of 4>j(r) made in Eq. 



Appendix B: Orthogonality Problem 

In this Appendix, we show that the eigenstates of the 
projected position operator X at finite size are not or- 
thogonal to each other due to the non-unitarity of the 
Berry connections. 



For simplicity, we focus on the case of a single occupied 
band; the problem is generic for any number of occupied 
bands. In this case, the Berry connection A x (k xi k y ) is 
Abelian, and the project position operator in the x di- 
rection is given by Eq. p5|): 



X — ^ ^ \k Xl ky) A x (k x •> k y )(k x -\- 1, ky 

I, (Bl) 



where the Berry connection is defined as [Eq. ( |26| )] 

A x {k x ,k y ) — ^ ^ 6 01 ^ x u rn a (k x , k y )u n ^ a (k x -\- 1, ky ) . 

a 

Generically, at finite size we have \A x (k x: k y )\ < 1; only 
in the thermodynamic limit does this become an equality. 

Since X is diagonal in k yi we can solve the eigenvalue 
problem at each k y independently. The N x eigenvalues 
of X at k q , are 



- i2KX/N "D(ky)X x (ky), 



(B2) 



labeled by X G [0 .. N x ). Here, the unitary phase factor 
^x(ky) is the N x -th root of the phase part of the Wilson 
loop along a fixed k y: defined in Eq. (30), while the posi- 
tive real number D(k y ) is the A^-th root of the absolute 
value of the Wilson loop, 



D(ky) 



N x 



J[ A X (ft, ky) 



l/N x 



(B3) 



The corresponding eigenstates are 



2itk x X / N x 

Tlkx A I 1 \ l^'^y/' 



\X,k y )=M kv J2e~ i 

k x 

and the normalization factor is given by 

[D{ky)] 2k * l 



-1/2 



(B4) 



Therefore, 



(X 1 ,k y \X 2 ,ky) = W ky ]*M ky 



y^ e t27rfc x (Xi- 



[D{ky)f 



Yl K x \A x (k, k y 



(B5) 



Evidently, in general the Wannier states localized in dif- 
ferent unit cells at each k y are not orthogonal at fi- 
nite size, preventing a nice mapping into the orthogonal 
LLL orbit als. The Wannier orbitals become orthogonal 
quickly in the thermodynamic limit due to the oscillatory 
nature of the sum, which averages to zero if X\ 7^ X%. 
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exist a Bravais lattice shift (A^, A^) G 1? and sublattice 

x y\ 



displacements (e§,e|) G M 2 such that 



- A 



a "I" fc /3 



Appendix C: Inversion Symmetry 

1. Inversion Operator 

We derive the inversion transformation of the expo- 
nentiated position operators x and y in Eq. (81 ) and pro- 
vide an explicit expression for the sewing matrix elements Note that (A^,A^) depends only on a, not on (3. Dif- 
ferent orbitals f3 may share the same value of (ej§, ejg). In 



-y-4 = -y-A»+6» 



(Cl) 



e <h x ,k y defined in Eq. ( [82] ). 

By properly choosing the reference point in each unit orbital space, inversion is implemented by 
cell, we can put the inversion centers on the Bravais lat- 
tice (x,y) G Z 2 . Inversion takes the a orbital located at 
(x + e%, y + e y a ) to (—x - e%, -y - e y a ). For the insula- 
tor to be inversion symmetric, we must have an orbital at 

(— x — e^, —y — ^a) before taking inversion. So there must where the matrix P is unitary and satisfies P 2 = 1. 



P\x,y,a) = 5> a/3 |-* - A*, -y - A^,/3), (C2) 

(3 



The position operator y transforms under inversion 



Y^ a p\-x - A*, -y - A», /?>(-* - -y - A^ )7 |P 7a 



EE' 



-i27v(-y-Al+el)/N y 



^2F a/3 \x,y,f3)(x,y^\. 



7Q! 



07 



J2J2e i27v(y+ ^ )/Ny 

x,y f3*y \ a 



\x,y,/3){x,y,j\ 



(C3) 
(C4) 

(C5) 
(C6) 



In the second ste p we changed the summation dummy index (x, y) —> {—x — A^, —y — A^), and in the third step we 
made use of Eq. (Cl). Similarly, we can show that VxV — x^. In momentum/orbital space, we have 



V\k x , k y ,a) = -±= e «*(W".+M/"») £ W af3 \-x - A£, -y - A» , /?) 

V^x^y x y p 



= e -i27r(k x Al/N x +k y A y a /N y ) ^ 



(C7) 
(C8) 



The sewing matrix elements 

e^--v=(k x ,k y \V\-k x ,-k y ) (C9) 
can be expressed in terms of P^, 

gi&a,,*,, = ^2^2 e i2ir(k x Al/N x +k v Ay/Ny) 

a (3 

x u*p(k x , ky)F al3 u a (-k x , -k y ). (CIO) 

2. Inversion Symmetry in Curvature Fluctuations 

We study the inversion property of the curvature fluc- 
tuation Uy(ky) = Uy(ky)/\Uy(ky)\, WllC^ Uy(ky) IS dC" 

fined in Eq. (53) as: 



N x 



Uy(ky) 



Wj^ky) 



1 



Nx^Wjk X ,ky) 



(53) 



First, we study the Wilson loop W^(k x , k y ). The inver- 
sion transformation of the Berry connection in Eqs. (86) 
and ( |87| ) relates W^(k x , k y ) to the Wilson loop around 
(N x — k Xj —k y )(N Xj —k y )(N Xj —ky — 1)(N X — k x ,—k y — 
l)(N x — k x ,—k y ), except for the corner case W_(0, k y ) = 
1. Due to the loop structure, all the sewing matrix ele- 
ments cancel completely, and we have 



W_(N X ,-k y ~l) 



1. 



(Cll) 



It is easy to check that the corner case k x = is incor- 
porated as well. 

Then we study the denominator W^(k x , k y ) = 
[i JL x{ky)) kx , with Hx{k y ) defined by Eq. (55). We can 
prove that 

W.(k X ,ky)=W_(k X ,-ky-l). (C12) 

Setting k x = N x in Eq. ( |C11| ), we have [recall that 
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W_(0, ~ky ~ 1) = 1] 

W_(N X ,ky)=W_(N X ,-ky-l). 



(C13) 



Using Eq. (|57J) to fix the branch choice, we can take the 
A^-th root of this equation and get 



l^xiky) — f^x{ ky !)■ 



(C14) 



Taki ng th e fc^-th power of the above equation proves 
Eq. JC12), 



Combine Eq. (Cll) with Eq. (C12), we can show that 
the curvature fluctuation U y {k y ) satisfy 



[Uy(ky 



= ly [WJk x ,k y )Y 



1 i \ 



W_(N X - k x 



1) 



WJN X , ~ky ~ I)/ W_(k X , ~ky ~ 1) 

WUN X -k X ,-ky-l) 
N X W.{N X) ~ky ~ \)/Wjk X , -ky ~ 1) 



l ■ — 

— T 



iv, ^ w_ 



W_(N X -k X ,-ky-l) 
W_(N X -k X ,-ky-l) 



— Uy( — ky — 1). 



(C15) 



Here we have used W^(N X , k y ) = W^(N X , k y ), which fol- 
lows from the definition of W_(k x ,k y ) = [/^ x (k y )] kx and 
Eq. ([56]) . This proves Eq. ([96]) ? namely, 



Uy(ky)U y (-k y ~ 1) = 1. 



3. Inversion Transformation of the Wannier States 

We prove the inversion transformation of \X,k y ) in 
Eq. ([99]). 

Using [X x (k y )] Nx = W x (k y ) [Eq. (30)], we can rewrite 



the inversion transformation of W x (ky) in Eq. ([88]) as 



i N x 



1. 



[K(ky)\ X (-ky)} 

Taking the A^-th root, we have 

X x (k y ) = e ^ 7rn / Nx j\ x ^— k y ). 



(C16) 



(C17) 



where n G Z is to be determined. Recall that by defini- 
tion, both X x (k y ) and X x (—k y ) must have their argument 
angles in (—2tt/N x , 0] mod 2ir. So there are the following 
two cases: 

"on-site" If W x (k y ) = 1, we have X x (k y ) — X x (—k y ) = 
1, and thus 



A^(^) = l/X x (-k y ). 



(C18) 



The Wannier centers are located on the Bravais lat- 
tice. 

"off-site" If W x (k y )^l, we have 



X x (k y ) — e * / x /X x ( k y ), 



(C19) 



since both X x (k y ) and X x (—k y ) have argument an- 
gle in (—2ir/N x ,0) mod 27r. In particular, for 
Wa^fej,) -1, we have A a; (fc 3/ ) = X x (-k y ) = 

„-ilv/N x 



From the inversion transformation of the Berry connection A x {k y ) in Eq. (86) we can get 



N x k x 



N x -1 



J^J A x (k, ky) I | (ft? ) — W A x (tv,k y ) — e ^ 



N x -k x ,k y g i£,N x ,k y 



k Ik K—N x —k x 

Therefore, we have the inversion transformation of the Wannier states 

N x rA /7 \ik x 



e^ve- l ^ x ^V\X,k y ) = -L V e- i2 * k * x ' N " ^ x{ky ^ - e ^,^-^,^ |_^ ; 



iV, 



2nk x X/N x i^x(ky)] 



N x k x 



Y\ K X x A x (k, k y ) 



k x 
N x 



e i«« x ,*„ -fc^.^y \k x , —ky) 



(C20) 



(C21) 
(C22) 



1 g c «x fc .x/AT. [A x (^)]-^ J 3:l^ K, f kv l. e^-^-^\k x ,-k v ) (C23) 



-i2irk x (-X)/N x i^x(k y )] 



' I ^CC 5 k y ) 



' e -i*y(-X,-k v ) \_X^ _ ky) if ^(^) = 1? 

e -i^(-x-i,-fcj |_ x _ 1? otherwise. 



(C24) 
(C25) 
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Here in the second step we performed a change of sum- 
mation variable k x —> N x — k x and made use of the in- 
variance of each term under k x — >• k x + N x , while in the 



last step we used Eqs. (C18) and (C19). 



4. Inversion Transformation of the ID Index 

We prove the inversion transformation rule of j x > k v in 
Eq. fll07|), namely, 



_jx,k y _ x if W^(0) = — 1 and N y odd; 



otherwise. 



(107) 



W e wil l use the inversion transformation of (X, k y ) in 
Eq. (106), namely, 



(-X,-ky) 

(-X-1,- 



ifW X (ky) = l, 

otherwise. 



A discussion of the location of W x (k y ) = 1 is thus neces- 
sary. At a generic k y , due to the absence of any protec- 
tion by symmetry, in general we can assume W x (k y ) ^ 1. 
With inversion symmetry, We can safely assume that 
W x (k y ) = 1 happens only at k y = or otherwise 
the flow of the Wannier center would exhibit "zigzag" 
patterns. 

We only elaborate the proof for C = +1; the situation 
when C = — 1 is completely analogous. We discuss the 
three cases of S y in Eq. (94) one by one, namely, 



if W x (0) 

N y /2 if W x (0) 

{(N y -l)/2 ifW x (0) 



-1 and N y even; 
-1 and N y odd. 

' = XNy + Cky + Sy IS 

Ny). As emphasized in 



Recall that the ID label j x > 
defined with Ck y + S y G [0 
Sec. |IIIC[ before performing the mapping, we need to 
first shift k y to the principal Brillouin zone (pBZ). The 

[0 .. N y ). The three cases 
the corresponding pBZ are 



pBZ is defined by Ck y + 6 y G 
for the shift parameter S y and 
illustrated in Fig. |5j 

In the first case, W x (0) = 1. In this case we have 5. 
0, and the principal Brillouin zone is pBZ 



k y = 0, we have 1^(0) = 1 and —k y 
thus under inversion 







[0 .. Ny). If 

G pBZ, and 



J 



XN V 



-XN V 



-J 



If k y G (0 .. Ny), we have W x (k y ) ^ 1 and N y -k y G pBZ, 
and thus under inversion 



XN, 



^(- X -l)Ny+Ny 



In the second case, 1^(0) 



have S y = and pBZ = [- 



Ny 

2 



-1 and N y is even. We 



AT, 



?-). Due to \C\ = 1, 



we have W x {-^fj = - 

we have W x (-^) = 1 and -N 



l/W x (0) = 1 [Eq. ([93])]. Ifky 

N v 



pBZ, and thus under inversion 

j x <~- =XN y 



2 j — x tAiivj. fry — 2 ^ 



2 2 



•>■-*>- 



If & y G (-^ .. we have W^/fy) 7^ 1 and € 

pBZ, and thus under inversion 



N„ 



(- X -l)Ny + (-ky)+^ 



~X,k v 



In the third case, W x (0) = — 1 and N y is odd. We have 



W x (k y ) ^ 1 and S y 



N v -1 



and pBZ 



N v -U 



(illustrated in Fig. [5^). For any k y G pBZ, we have —k y G 
pBZ, and thus 



XNy + ky 



N, 



(- X -l)Ny + (-ky) + 



N v 



-J 



X,k v 



1. 



Su mma rizing the three cases, we have proved 
Eq. (107). We emphasize that for a given system [and 
thus given W x (0) and N y ], all the j x > k y indices obey 
the same transformation rule, but the trans forma tion of 
(X, k y ) indices still comes in two cases [Eq. (|106|)]. 



5. Inversion Symmetry of the Many-Body States 

We study the inversion operation on the many-body 
states s,r)iat- From Eq. pg} we have 



({X,k y }r\^;s,r) lat = 

{{X, k y }\V\{X u ky tl }) ({j x ^m s, r). (C26) 

Plugging into the above equation the inve rsion transfor- 
mation of the Wannier states in Eq. (105), the sum over 
(Xi, fc^i) reduces to a single term, with (Xi, k y ^) equals 
either (—X, — k y ) o r (—X—l,—k y ), depending on the con- 
dition in Eq. ( |105| ), namely, whether W x (0) = — 1 and N y 
is odd. For both cases, the Wannier-basis inversion ma- 
trix element has the same value 

{{X,k v }\V\{X u k v +}) = e^°[\ y (0)u v ]- 2k y. (C27) 

And we want to relate the FQH amplitude factor 
({j x ^ k y^}\^f;s,r) back to ({j x > ky }|^; s, r), but this de- 
pends on the inversion transformation j x > k v — ». j x i, k y,i m 
As shown in the previous section of t his Appendix, we 
need to discuss two cases in Eq. (107) for the inversion 
transformation of the ID index j x ^y . 



29 



First, consider the case where W x (0) = 1 or N y is 
even. In this case, j x > k v — >. j x i> fc y,i = —j x ^ k v under in- 



version. Also, according to Eq. (94), the shift 5 y is either 
or Ny/2 , and thus N y divides 2S y . From Eqs. (fC26) 



and (C27) we then have the amplitudes of the inverted 
state 

({X,k y }\V\^s,r) lat (C28) 
= e iN ^ [Xy(0) uj y ]~ 2 E *v <{-j*.*v}|* ; s , r ). 

The inversion transformation of the FQH state in 



Eq. (80) then relates this to the amplitude of another 



FQH state s,r), and we can map it back to the cor- 
responding state on lattice f)i at : 



({X,fcJ|P|M/; S) r) lat 



(C29) 



= C* e l ^«°-° [A„(0) oj y }- 2 ^ fe » ({X, fc y }|M/; s, f) lat . 



The indices s, f are defined uniquely by Eqs. (78) 
and fl7\ . The constant is the eigenvalue correspond- 
ing to~|\fr) of a composite operator [Eq. ( [75] )] and it takes 
value from ±1. The inversion operation thus relates the 
components of the many-body state s, r)i at with total 
momentum [Eq. (69)] 



(K x ,K y ) = [K x + sN e ,CK y + C(r-5 y )N e ] mod (N x ,N y ), 
to those of the state s, f)i at with total momentum 



(i^, K y ) = [k x + 57V e , Ck v + C(f - 5 y )JVc 
= (-K x ,-K y ) mod (N x ,N y ). 



(C30) 
(C31) 



thanks to Eqs. ([78]) and ([79]) and the fact that N x divides 
Nq x and N y divides 25 y as noted earlier. 

The coefficients of the inversion transformation in 



Eq. (C29) appear to vary for each component, as ky 
is different for different components due to umklapp 
(mod N y ) processes. However, since with inversion 
[X y (0) uj] 2Ny = 1, according to Eq. ( [98] ) ? the troublesome 
factor [X y (0) uj y ] 2 ^ k y depends only on Yl k y m °d N y , 
and thus all the non-zero components ({X, k y }\^;s, f)i at 
have the same value of Y ky mod N y , na mely, — K y . This 



means that the coefficients in Eq. ( C29| ) are actually the 
same for each component. Therefore, the inversion op- 
erator V indeed brings one state to another within the 
g-fold multiplet, 

P|*;s,r)iat (C32) 
= C^e^^'°[A,(0)^] 2C ^+^-^)^ |*;5,f>ut- 

We now proceed to the other case where W x (0) = — 1 
and N y is odd. The majority part of the above argu- 
ment carries over, except that now j x > k v is mapped to 
—j x ' k y — 1 under inversion. The —1 can be pulled out 
as a translation Ti*L and absorbed into the state s, f) 



In the current case we have S y = —1 — S y mod N y . and 
thus the state s, f — l)i at has the total momentum of 
the inversion partner of s, r)i at , 

(K x , K y ) = [k x + sN e , C Ky + C(f - 1 - S y )N e ] (C34) 
= [-K x -sN e ,-C Ky -C(r-5 y )N e ] (C35) 
= (-K x , -K y ) mod (N x , N y ). (C36) 

Again, the inversion operator V brings one state to an- 
other within the q-fold multiplet, 

P|*;*,r)iat (C37) 
= ^e iN ^"[Xy(0)LO y } 2C ^ r - S y^ |*;5,f- l) ut . 

The index f — 1 needs careful handling, since |\£;s, f) 
acquires a phase whe n f - > f+q/ q x , as shown in Eq. (19). 

C ombi ning Eqs. (C32) and (C37), we have proved 
Eq. flTM]). 



Appendix D: Recombination of the LLL Orbitals 

In this Appendix we prove the recombination formula 
for the LLL orbitals: 



oH 4 2 / 



i2-Kjm/N ( f ) 



(Dl) 



where the orbital wave func tion (/) m (x ,y) i s given by 
Eq. (|3|), and by Eqs. ( [Tig ) and ( [118] ). 

Consider the single-particle magnetic translation oper- 
ator T{Li/N(h) in the Landau gauge A = Bxe y . This 
operator commutes with T(Li) and T(L 2 ), and thus is 
compatible with the periodic boundary conditions in the 
LLL. From the explicit expressions, we can find the ac- 
tion on the Landau orbitals 

TiU/N^jfay) = <f> j+ i(x,y), (D2) 
TiU/N^ix, y) = e i2 ^/ N ^ V j(x, y). (D3) 

The second equation shows that the operator T(Li/A^) 
has N^ distinct eigenvalues within the LLL. Since the 
LLL as a Hilbert space has dimension TV^, each of the 
eigenspaces of T(Li/iV^) within t he L LL has to be non- 
degenerate. Moreover, from Eq. (D2) we can construct 



another state that has the same eigenvalue as <j) v Ax,y), 



T{U/N^Y,^ i2 ^ /N ^rn{x,y) 



_ ^nj/Nt y^^ p -i27rjm/N (f>f 



>&y)- (D4) 



see Eq. (JTSJ)] . The amplitudes of the inverted state reads Therefore, we must have 



{{X,k y }\V\*;s,r) lat (C33) 
= C* e iN ^°-° [X y (0) La y ] " 2 £ k * ({X, ky}\9;8,f- l) lat . 



b v j(x, y)(xJ2 e- i2 ^ jm/N ^ m (x,y). (D5) 
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The proportionality constant cannot depend on (x, y). It 
can be fixed by comparing the two sides at (x, y) = (0, 0) . 
We find 



1. C7 = +l 



$(0,0) 
0m (0,0) 



1 



1 1 



(D6) 



Comparing Eqs. (69) and (122), we find that the state 
l^sVOiat with ( s Vn given by 



m I i 



), (D7) 



s' = s - S x , 



r = r — S u 



(El) 



has the same total momentum a s the state s, r)i at . 

/-z — j(ir/A-f)/ t >\ i o / i \ • n ^i-i The above equation replaces Eq. (125) in the main text. 

where £ = ^LxjL^ e H7r/4 and #3(z|t) is the third H ^ m v| p 

We now examine the corresponding FQH states. Sub- 
stituting in Eq. (17) gives 

(D8) \y;s',r') = e^(Tyj s *(T c x J- s y\y;s,r), (E2) 



Jacobi-$ function, 



M*\t) = ^ e <W+i27rn *, Imr > 0. 



Using the discrete Fourier transform formula in Ref. [551 where e ^ = e i27rS ^ +rNe ^ N ^ is an inconsequential 



find 



N 6 



phase factor. We find the ac tion of T£ m and T^ m on \{j}) v 
by gauge transforming Eq. (115): 



i _ 

? 9„ C J-\iSl) — - -il-Kim/N^ a i ml t \ 



(E3) 

The proportionality constant at (x, = (0, 0) is then Therefore, 
found to be e l ^l^~ e l 2 ^ / y^N^. This pro ves the recombi- 
nation formula of the LLL orbitals Eq. ( |Dl| ). v ({j}\^; s ,r) = e^e~ i2n ^ j ~ 6 ^ 5y/N<i> v ({j-S x }\^-,s,r). 

' (E4) 

We now examine the lattice amplitudes 
({k x ,k y }\i&is, r)' obtain ed by the Wannier mapping. As 
explained in Sec. VII C we put all the momenta {k x , k y ) 
in the principal Brillouin zone determined by 



Appendix E: Proof of Quasi-Isotropy 



In this Appendix we generalize the proof in Sec. |VII C 



to cover C = ±1 with generic values of S x and S y . The 
only hurdle compared with the solved case with C = 1 
and S x = S y = is that the states s, r)i at and 
g, r)[ t reside in different momentum sectors [Eqs. ( [69] ) 
and (pgf]. 



-fex + <k G [0 .. AT,), ky -\- Sy E [0 .. Ay. (E5) 



We discuss the case of C = 1 and C = — 1 separately. directly. 

I 

We find 



This enables us to use the Wannier mappings j x > k v and 



v ({tf k *}\*;s',r') = V ({YN X - k x + 5 x }\*;s',r') = e ip e- i2 ^ YN ^ k ^/ N \{{YN x - k x }\*;s,r). (E6) 
Plugging in the Fourier transform in Eq. ( |12Q[ ) gives 



!__ ^ e~ i27T E ^ ei2n E ZyY/iVy e_i27r E fea: '* /JV * ({ }|#; s, r) . (E7) 



{*} 



Here we have rewritten each sum over m in the Fourier transform as a double sum over all unit cells X and all 



l y points in the principal Brillouin zone, and e v 
({k x ,k y }\^;s f ,r f )[ at becomes 



-Wc(f-f) 



2 Ms a constant phase factor. Then, the core part of 



1 N ^2c~ i27rT ' KY/Nv v ({jv ,kx }\^]s\r f ) = c il3+i ~ f ~ i2 ' K Tj kx K /Ns 1 N e ~ i2ir ^ k * x l N * ({j x > k y }|^; s, r) . 



y {y} 



Plugging the above equation into Eq. (123) gives Eq. ( 128) with an inconsequential extra factor e 1 ^ . 
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2. C 



Therefore, the FQH amplitudes of the primed states are 
given by 



In this case, we need to find (V, r') such that 

k x + (V + S x )N e = -k x - sN e mod N x , 

K y + r'N e = — (Ky + (r — 5 y )N e ) mod AT^. 



(E9) 



Putting in Eqs. ([78]) and Eq. ([79]) gives 
s' s-S x y = r + 5 y , 
where (s, f) are related to (s, r) by 

({-i}|*;s,r) = C*<{j}|*;«,f>. 



(E10) «({?}!*; *'y> = e^e- <2 '2:«-' 5 -) a »/ Jv ^{i-<J x }|*; a ,r>. 

(Ell) 

We emphasize that the above steps requires the inversion 
symmetry of the FQH states only; there is no requirement 



(80) on the FCI system. 



J 



The following development largely parallels the C = +1 case. Putting all (k Xl k y ) in the principal Brillouin zone, 
we find 



v ({tf k *ms',r') = Q l e^e +i2 ^ YN -- Ck ^/ N \{{-YN x + Ck x }\^;s,r). 
Plugging in the Fourier transform in Eq. ( |12Q[ ) gives for C = — 1 



(E12) 



v ({fr k °ms , y) = Q 1 e 



1 i/3+ij _ 



1__ J2 e i2nC E k * x / N * e ~ i2nC E ^ y/A V 27r E k * l v/ N * ({j x ^ }|#; 5, r). (E13) 



{X} 



{ly} 



And the core part of ({k x , k y }\^; s' , r')| at becomes 



/TV, 



e ~ i27V ^ kyY / Ny v ({j Y,kx }\ty; s', r') = ^— e - i2 ^ C E Kky/N^ ^ e ~ i27r E k x X/N x ({j X ' fc ?/}|^r; s? r ). 

(E14) 



Pl uggin g the above equation into Eq. ( |123[) gives 
Eq. (128) with an inconsequential extra factor e 2 ^/(^. 



Appendix F: Continuum Limit 

We consider the continuum limit of a Chern insulator, 
N x — >• 00, A^ — >• 00 with a fixed aspect ratio N x /N y . 
For simplicity, in the following we focus on the case of a 
square Bravais lattice, with bi • b2 = and |bi| = |t>2 1 - 
The formulas in the generic case can be obtained without 
much extra effort. We define the continuous coordinates 
x G [0,N X ), y G [0,N y ), and the continuous momenta 
k x ,k y G [0,2tt), 



k x = 2irk x /N x , 



k y = 2irk y /N y . (Fl) 



Normally, the continuum limit involves sending the unit 
cell size to zero. To make direct comparison with the 
formulas in Ref. [I5j however, we do not take the limit of 
the unit cell size going zero, and work instead in a unit 
system in which |bi| = 1 1>2 1 = 1- 



1. Berry Connection 

We show that the discrete connections defined in 



Eq. ( 26 ) reduce to the Berry connection in the usual def- 
inition in the continuum limit. 

Recall that the Bloch state is given by 



\k x ^k y ) — ^ ^ u a (k xi ky)\k x , ky , ol) 



(F2) 



= ^_ E e^ x+i ~ k ^u a (k x , ky)\x, y, a). 

V x y x,y a 

(F3) 

The single-particle orbitals are defined by the real-space 
wave functions. 



(x, y\x, y, a) = (j> a {x - x - e%, y - y - e y a ). 



(F4) 



Here the function (f) a (x, y) is centered around x = y = 
and we have made explicit the sublattice shift e% y . 

We assume that that the orbitals are orthonormal in 
the sense that 



(x,y,a\x',y',P) = S a( 35 xx >5 yy > . 



(F5) 
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And in the spirit of tight-binding approximation, we fur- 
ther assume that the orbitals are localized enough such 
that 



(x,y,ot\x\x',y',a) = (x + e*)5 a p5 xx >5 yy > , 
(x,y,a\y\x',y',ot) = (y + e y a )5 al3 5 xx >5 yy > . 



(F6) 



We can calculate the Berry connection a = (a x ,a y ) 
using the usual definition 

&(k x ,k y ) = -i(u %a $jV % \u %w $ y ), (F7) 

where the periodic part of the Bloch state is given by 

K x ^=e-^-^%A>- (F8) 

The result is 

a x(kx-> k y ) = — ^ ^ u a (k x , ky)(idfc + e^)u a (fc £C , ^^), 

a. 

a y{k X i k y ) = ] u a(^xi k y )(id^ ^a) u Oi(]^x^ ky) . 

a 

(F9) 

We assume that the Bloch state \k x ,k y ) is periodic in 
the fca direction, and satisfies 



(F10) 



in the k y direction (similar to the Landau gauge). The 
Stokes' theorem is applicable, but the non-zero Chern 
number does not obstruct the construction of a smooth 
gauge, as we have abandoned the periodic gauge condi- 
tion in the k y direction. This is exactly the same bound- 
ary condition adopted in Ref. [15j 

Since the exponentiated Berry connection restores uni- 
tarity in the continuum limit, we can relate its phase to 
a.(k x ,k y ). To order 0(1/N X ), we have 

A x {k Xl k y ) — ^ ^ 6 a ^ x u a (k Xl k y )u cx {k x -\- -j^— , k y ) 

a 

~ x ~ ^Yl^^xJ^y)^^ +e x a )u a (k x ,k y ) 

x a 

^ e i2ira x (k x ,k y )/N x ^ (Fll) 

Similarly, to order 0(l/N y ), 

Ay(k X ,ky) « e ^ y (k X ,ky)/N y ^ ^ 12 ) 

In the thermodynamic limit, the above approximate re- 
lations become exact. 

2. Wilson Loops 

Define the phase angle (polarization) 

0(k y ) = / dp x a x (p x ,k y ) (F13) 
J o 



The Wilson loop along a fixed k y is given by 

W x (k y ) = e i0 &y\ (F14) 

The Bloch state \k Xl k y + 2tt) is connected to the state 
\k x ,k y ) by a gauge transform. Specifically, since we as- 
sume a y (k Xl k y ) = 0, we have 

— exp z / dj9y {kx i Py ) (^cc i ky -\- 2tt | k x , /c^ ) 
= (k x ,]^ + 2v\k x ,ky). (F15) 
The gauge transform is thus given explicitly by 

\k x ,k y + 27T) = [Wy(k X )]*\k X ,ky). (F16) 

The Wilson loop W x (k y ) is gauge invariant. Naively, 
one might expect that the phase angle 0(k y ) is also peri- 
odic. However, in t he con tinuum limit, this is not true. 
As explained in Sec 



IIIC 



when k y gradually increases by 
27r, W x (k y ) winds around the unit circle at the origin of 
the complex plane. The winding number in the clockwise 
direction is given by the Chern number C. Since we as- 
sume the smooth gauge for the Bloch states, 9{k y ) must 
be smooth as well. Therefore, we have the quasi-periodic 
condition^ 



6(k y + 2tt) = 0{k y ) - 2nC. 



(F17) 



In the light of Eq. (F16), this can also be understood as 



(F18) 



(F19) 



6(k y ) 6(k y ) - 2ttC 
under the gauge transformation 

\k x ,k y ) [W y {k x )Y I k x 5 ky) . 

3. Maximally Localized Wannier States 



In the continuum limit (which unfortunately cannot be 
used to compute overlaps with the exact ground states 
from numerical diagonalization on the finite-size lattice) , 
the projected position operator X = PxP becomes uni- 
tary. Therefore, the eigenstates \X,k y ) of X (the unitary 
part of X) are just the eigenstates of A', the ID maxi- 
mally localized Wannier states. The change in the gauge 
condition in the k y direction does not affect the form of 
the Wannier states, as the construction is performed in- 
dependently at each k y . We can rewrite the definition in 
Eq. (32) using the continuum variables, 
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\X,ky) 



/N x ± 

k. 



^exp 



ik x X + ik x - i / dp x a x 

In J 



(F20) 



r 



This expression is exactly the same as in Ref. 15. The 
center position of this state can be related to the phase 
angle by 0{k y ) 



(X,k y \x\X,k y )=X 



g(fcy) 
2tt 



(F21) 



This is fundamentally different from the transformation 
in the periodic, non-smooth gauge on a lattice [Eq. ([39])] : 
in addition to acquiring a phase, the Wannier state may 
hop to a different unit cell upon a gaug e transformation. 
The transformation given in Eq. ( F26 ) differs from that 
of Ref. HSl 



We now move on to discuss the gauge transform of 
these Wannier states. This is where the current discus- 
sion diverges from Sec. [ITT] and Ref. 15. Consider a gauge 
transformation 



The discrepancy between Eq. (139]) and Eq. ( F26 ) comes 



from the smooth gauge condition. Due to the winding 
motion, the argument angle of W x (k y ) is multi- valued by 
nature. At finite size, we require X x (k y ) to be the A^-th 
root of W x (k y ) with argument angle in (—2tt/N x , 0]. This 
corresponds to confining the argument angle of W x {k y ) 
to the branch (— 27r,0]. In the continuum limit, 



\k x ,ky) -> e^ k ^\k x ,k y ). (F22) 
The Berry connection transforms by 

) + V%v(k x ,ky), (F23) 

and thus 

/ dp x a x (p x ,k y ) / dp x a x (p x ,k y ) 
Jo Jo 

+ V(k x ,ky)-v{0,ky)' (F24) 

In particular, for k x = 2tt, 

6{k y ) — >> 0(k y ) + t](2tt, k y ) — 77(0, k y ). (F25) we use the periodic part of |X, k y ) and calculate 



[K(k y )] kx = [W x (k y )] k ^ N * -+ [ e ^(^)]^/( 2 ^. (F27) 

The smooth gauge condition in the continuum requires 
0(k y ) to be smooth as well, and thus makes it impossible 
to fix 0(k y ) to (-2tt,0]. 

4. Curvature Fluctuations 



We now derive the connection between Wannier states, 



the continuum analogue of Eq. (52). Similar to Eq. (F7), 



This is consistent with Eq. (F18). Putting all of these 
together, we find that 

\X,k v ) e^°^\X- ±[ V (2n,k y ) - v(0,ky)], k y )- 

(F26) 



a v (X,k y ) = -i(X,k v \e ik y y d~ k e- ik y y \X,k v ). (F28) 



From Eq. (F20), we obtain 



a y (X,k y ) =d ly $ y (X,k y ) + — 



dk r 



dz 0{k x , k y 



2 - , fe y v v"-*< "•»' / dp x a x {Pxi ky) + a,y(k x , k y ) 



(F29) 



Here we have used the limit Y^k* ~^ 2^ Jo* ^ x ^° ^ e enc ^ U P 
make the notations coherent. Recall that a y (k xi k y ) =0, 
the third term in the bracket vanishes, while the first two 
terms can be related to the Berry curvature 



fxy{k X jk y ) — dj: a x {k x ^ky) d^, a y (k x ,k y ). (F30) 



a v (X, k y ) = d iy $y(X, k y ) 
+ — dk x dp x f xy (Px,ky) - f xy (k y ) , (F31) 



2tt 
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where the average Berry curvature f X y(ky) is defined by 



- ~ _ i r 2 * ~ ~ ~ 

fxy(ky) = J ^k x f xy (k Xl 



)• (F32) 



Similar to the discrete case, the fluctuation of Berry cur- 
vature appears in the connection between adjacent Wan- 
nier states, even more sharply than the finite-size formula 
in terms of W^(k x , k y ) and W_(/%, k y ). 

We stop the exposition of the continuum case here, as 
the rest of the construction could be obtained directly by 
taking the continuum limit of the finite-size formulas. 



Appendix G: Comments on Qi's Prescription 

In this Appendix we discuss issues in Ref. [15] that we 
have found while trying to extend the prescription to the 
torus geometry. 

The original proposal was to map the eigenstates of 
the projected position operator in the gauge a y = to 
the states in the LLL, and use these Wannier states to 
build the cylinder FQH states on a lattice. One would 
expect a similar construction would hold for the torus 
FQH states. As pointed out in the main text, the subtle 
points are that we need to apply the phase fix discussed in 
Sec. |III A] to restore orthogonality of the Wannier states, 
and that we have to modify the Wannier-LLL index map- 
ping according to Eq. (43). 



Based on the observation that the Wannier center shift 
by one unit cell when k y increases by 27r, Eq. (3) in 
Ref. H5l claimed that \X,k v + 2ir) = \X + C,k y ). However, 
as pointed out in Appx. |F 2\ we can t hink of k y —> k y +27r 
as the gauge transform in Eq. (F19), and thus the Wan- 
nier state \X, k y + 2tt) is given by Eq. {F26), 



|X, ky + 27T) = [Wy(0)]* \X + C, ky). (Gl) 

Ref. H51 missed the extra factor [W y (0)]*. 

This seemingly innocent problem is actually catas- 
trophic, unless W y (0) = 1. It is directly responsible 
for the low overlap between the constructed many-body 
state and the ground st ates f or several important FCI 
models examined in Sec. 



VIA 



Assuming \X,k y + 2tt) 



\X + C,k y ) as in Ref. [I5j we can limit k y to a single 
Brillouin zone when building the N x N y Wannier states. 

Consider the simplest case with Chern number C = +1 
and the shift parameter S y = 0. The principal Brillouin 
zone is given by (k x ,k y ) G [0 .. N x ) x [0 .. N y ). We 
need to implement of the gauge condition a y = on the 
lattice. The lattice analogue of a y = is setting the 
exponentiated connection A y = 1. As pointed out in 
Sec. |IIIB[ we only need to enforce this condition along 
k x = 0. Namely, we should fix the phase of |0, k y + 1) 
relative to \0,k y ) by requiring A y (0,k y ) = 1. Starting 
from the Bloch state at (k xi k y ) = (0,0), this can be 
done recursively in the interior of the principal Brillouin 
zone till we reach the boundary (k xi k y ) = (0,N y — 1). 
This introduces a jump in the phase of A y (0, k y ) at the 



boundary of Brillouin zone. The Berry connection along 
k x = is given by ^(0, k y ) = A y (0, k y + N y ) and 



A y (0, k y ) — 



1 if ky G [0 .. Ny-1), 

Wy(0) liky = Ny-l. 



(G2) 



The value of A y (0,N y — 1) reflects the relative phase of 
the Bloch state at (0, N y ) as fixed by the periodic bound- 
ary and the (gauge- invariant) Wilson loop. We empha- 
size that A y (0,N y — 1) is the connection between the 
Bloch states at momenta (0,N y — 1) and (0,0), rather 
than the states at (0, N y — 1) and (0,N y ); the value of 
A y (0, N y — 1) in Eq. (G2) is not in conflict with a y = 0. 
In this particular gauge of Bloch states, the original 



construction amounts to Eq. (32) with e l ® y ( x,ky ^ = 1. 



The analysis in Sec. IV C shows that for this choice 
of e l ® y ( x,ky \ the resulting many-body wave function 
is indeed translationally invariant, in agreement with 
the claim in the original paper. 15 However, the connec- 
tion between adjacent Wannier states (X, k y \Y\X' , k y ) oc 
A y (0, k y ) changes abruptly at the Brillouin zone bound- 
ary k y = N y — 1. Unless W y (0) = 1, this strong k y de- 
pendence sets the Wannier states apart from the Landau 
orbit als, as shown in Sec. |IVB[ and destroys the FQH 
character of the constructed lattice trial state, in accor- 
dance with the numerical results presented in Sec. VIA 
We now check the inversion symmetry. Putting 

e i<S>y(X,ky) 



1 in Eq. (99), we have 



e*°' k vP\X, k y ) = 



-X — l,—k y ) otherwise. 



(G3) 

Recall that the sewing matrix elements e z ^°> ky are re- 
lated to the Berry connections. Plugging Eq. (G2) into 
Eq. d87b, we find 



p*£o,o 



if k,, = mod N, 



Wy(0)e^ ' otherwise. 



(G4) 



Restricted by the inversion symmetry, W y (0) can be ±1. 
If W y (0) = 1, then all the analysis of our proposal in 
the main text can be reused to show that the original 
prescription in this case preserves the inversion symme- 
try on the many-body level. If W y (0) = —1, however, 
the constructed many-body wave function breaks inver- 
sion symmetry. The reason is very simple: since \X,k y ) 
receives a different coefficients upon inversion depend- 
ing on whether k y = mod N yi different components 
of the many-body wave function transforms differently 
under inversion. We expect the degree of the inversion 
breaking to decrease as the system size increases. 

It should be obvious that the above problems in origi- 
nal proposal is not limited to the case C = 1 and S y = 0; 
it erupts whenever W y (0) = — 1. 

Finally, we note that in a related paper, 22 a slightly 
modified formula for the ID localized Wannier states 
were proposed for a generic gauge condition on a y . In 
our notations, the updated formula reads 



35 



\X, k y ) = exp 



ik. 



v 2tt 



JO 



dpy a y (0,p y ) 



exp 



ik x X + ik, 



0(ky) 



JO 



dp x a x (p x ,k y ) \k x ,k y ). (G5) 



Here the real number 6 y is defined by 

r 27T 



'v - I dp y a y (0,p y ). 
Jo 



(G6) 



The exponential prefactor in Eq. ( G5 ) was introduced to 



fix the gauge choice between the Wannier states at differ- 
ent k y . 22 Judging by appearance, the formula may look 
similar to ours [Eqs. ( [32] ) and Eq. (63)] at first sight. 
However, further scrutiny reveals a fundamental differ- 
ence. The continuum formalism in Refs. [15] and [22] re- 
quires a smooth gauge that is periodic in the k x direction. 
Stokes' theorem then dictates that the gauge cannot be 
periodic in the k y direction. In particular, a y (0,p y ) is 
not periodic. Therefore, the phase e z0y is not the Wil- 
son loop along k x = and is not gauge invariant. We 
can actu ally remove the added exponential prefactor in 
Eq. (G5) altogether, by going back to the a y = gauge. 



Therefore, the updated formula for the Wannier state in 
Ref. [22] is actually no different from the original one in 
Ref. [I5j and they suffer from the same problem detailed 
earlier. 



Appendix H: Overlap between Random Vectors 

Consider two random unit vectors sampled indepen- 
dently and uniformly from the vector space C d . We are 
interested in the probability distribution of the absolute 



square of the inner product (the overlap) between them. 
We seek an expression for the probability Pd(u) for the 
overlap to be larger than u. 

Since the absolute square does not see the relative 
phase between the two vectors, we only need to study 
R d . The uniform sampling of unit vectors in C d corre- 
sponds to the uniform sampling by area over S d , the unit 
(d — l)-dimensional sphere surface embedded in We 
build a Cartesian coordinate system #2? • • •> x d) with 
the x\ axis pointing in the direction of one of the two 
vectors. Then we can see clearly that Pd(u) is just two 
times the surface area of the hyperspherical cap on S d of 
height (1 — y/u), divided by the surface area of S d . Using 
the hyperspherical polar coordinates, we find that 



Pd(u) 



2S d - 



i / dx(l-x 2 ^ d - 3 ^ 2 , (HI) 



where S n is the surface area of S n , given by 

i/2 



r(f + i) 



(H2) 



This can be evaluated using numerical quadrature. We 
note that at large d, the probability Pd{u) to have overlap 
higher than u is exponentially small, bounded from above 
by 



0[(l-u) 



(d-l)/2i 



(H3) 
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It is always possible to line up L2 with the e y axis for a 
torus of any aspect ratio. 

Notice that the original treatment in Ref . [16] put an ad- 
ditional factor of ( — l)^(^ e -i) m the formula, but this is 
not necessary for our purposes: (K x ,K y ) are defined by 
Eq. Q. Our choice can be regarded as an alternative la- 
beling aimed to simplify the formulas in this paper, at the 
cost of the identity of (k x , Hiy ) as true momentum. 
We emphasize that the operators S x and R y are defined 
in the continuum. In contrast to the discussions in Ref. 9, 
here we do not add a lattice pinning potential to the FQH 
setup. Were we to do that, S x and R y would not be le- 
gitimate translation operators as they contain one-body 
translations that are a fraction of a unit cell size. We bridge 
the FQH and the FCI sides through the properties of the 
second-quantized amplitudes, without directly migrating 
the continuum translation operators and their algebra to 
the lattice. 

A notable special case is when q divides N y . Since 
GCD(N y ,q) = GG~D(pN y ,q) = N y N 0x /N = q/q x , in this 
case we have q x = 1, and thus S x = (T^)* and R y — T c y m . 
Hence, the original g-fold states diagonal in (T^ el ,T^ m ) are 
already the simultaneous eigenstates of S x and R y without 
any recombination. 

We thank F.D.M. Haldane for stressing the importance of 
this point to us. 

The non-unitary connection A a {k yi k y ) (a = x,y) may 
vanish on a very small lattice. This problem is most pro- 
nounced on an N x x N y — 2x2 lattice with inversion 
symmetry. In this case, all the four Bloch states are in- 
version eigenstates. The connection A a between two Bloch 
states vanishes when the two states belong to different in- 
version eigenvalues. We do not consider such pathological 
cases in our treatment. 

Arguably we could also use the aspect ratio as a variational 
parameter. We leave this for future work. 
The absolute value of (X, k y \Y \X' , k' y ) does not depend on 
the choice of the phase e ^ (x ' fey) . It is given by \ U y {k y )\. 
In the case of many occupied bands, the inversion oper- 
ator is a matrix, but in the absence of additional sym- 
metries in 2D, it can always be represented in a diagonal 
[U(l)} Occupied formP 

The two-orbital model is not included in the 6x3 calcu- 
lations, because this model does not have a gapped topo- 
logical ground state at this lattice size. 
Notice that this quantity cannot be generalized to the case 
with more than one state in a momentum sector, for ex- 
ample the Laughlin FCI state on a 6 x 3 lattice, as the 
absolute value operation breaks the orthogonality between 
the degenerate states at the same momentum. 
Even when the C4 symmetry is present in the hopping 
amplitudes locally, the lattice aspect ratio has to be fixed 
to specific values (square, for example) for the symmetry 
to survive at finite size. 



